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I .  INTRODUCTION 

For  the  past  several  years  exploratory  development  bottom  interaction 
studies  at  ARL:UT  have  concentrated  on  the  acoustic  mechanisms  of  bottom 
interaction.  Substantial  effort  has  been  directed  toward  the  study  of 
reflection,  refraction,  and  absorption  within  realistic  models  of  marine 
sediments.  These  studies  have  had  as  their  principal  technical  goals  the 
identification  of  all  processes  important  to  low  frequency  acoustic 
propagation,  and  the  quantification  of  the  relative  importance  of  these 
processes.  These  studies  included  theoretical  investigations  complemented 
by  close  interaction  with  data  analyses  carried  out  in  concurrent  advanced 
development  programs  at  ARL:UT. 

During  the  past  contract  year  (FY79)  there  have  been  four  lines  of 
investigation:  (1)  lateral  subbottom  variations,  (2)  sloping  bottom 

effects,  (3)  subbottom  roughness  scattering  loss,  and  (4)  bottom  inter¬ 
action  sensitivity.  The  first  two  topics  have  been  grouped  together  for 
reporting  purposes  under  the  heading  of  range  variable  environment  studies. 
Of  these  four  study  areas,  the  first  two  are  expected  to  continue  into  the 
next  year,  whereas  the  second  two  are  terminated  with  the  work  reported 
here. 


This  year's  research  has  resulted  in  extensive  documentation, 
primarily  in  the  form  of  journal  articles,  some  of  which  have  only 
recently  been  submitted  for  publication.  Because  of  this  extensive 
reporting  the  present  report  contains  only  a  summary  of  the  principal 
results  together  with  references  to  the  papers  and  reports  where  complete 
documentation  is  to  be  found.  Three  of  the  most  recent  papers  are 
included  as  appendices  to  this  report  since  their  appearance  in  the 
literature  will  occur  long  after  this  contract  year  is  terminated.  In 
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addition  to  summarizing  the  work  carried  out  during  the  past  year,  this 
report  also  contains  bibliographies  of  the  complete  documentation  for  the 
contract  year  as  well  as  complete  project  documentation  to  date. 


II.  PROPAGATION  IN  A  RANGE  VARIABLE  ENVIRONMENT 


The  work  concerning  propagation  in  a  range  dependent  environment  is 
summarized  in  this  section.  This  work  is  subdivided  into  sloping  bottom 
and  subbottom  lateral  variability  problem  areas. 

\ 

The  primary  goal  of  the  work  on  lateral  variability  has  been  to 
determine  whether  lateral  subbottom  variability  that  might  reasonably  be 
expected  to  occur  in  marine  sediments  is  a  significant  factor  influencing 
acoustic  propagation.  In  long  range  acoustic  propagation  some  degree  of 
lateral  subbottom  variability  is  always  present.  The  sources  of  this 
variability  are  differential  sediment  depositional  processes  and/or 
propagation  over  differing  geophysical  provinces.  The  prime  focus  of  the 
the  lateral  variability  work  is  to  determine  whether  range  variability  of 
the  bottom  is  likely  to  be  significant,  and  if  so,  how  best  to  include 
it  in  the  description  of  acoustic  propagation,  possibly  through  average 
bottom  loss  or  average  geoacoustic  model  descriptions. 

The  main  goal  of  the  work  on  sloping  bottoms  is  to  determine  how  best 
to  describe  acoustic  propagation  in  these  areas  and  what  level  of  detail 
is  required  to  adequately  characterize  the  acoustic  propagation.  Sloping 
bottom  problems  of  particular  interest  concern  propagation  over  continental 
slopes  and  seamounts.  Another  goal  is  to  determine  the  role  of  the  bottom 
properties  of  the  slope.  In  other  words,  are  the  geometrical  properties 
of  the  sloping  bottom  the  driving  mechanism  influencing  sound  propagation, 
or  do  the  geoacoustic  properties  of  the  bottom  play  an  important  role? 
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A. 


Review  of  Previous  Work 


1 .  Lateral  Variability 

At  the  beginning  of  this  contract  year  the  work  on  lateral 
subbottom  variability  had  progressed  to  the  point  of  being  able  to 
estimate  the  importance  of  mode-mode  coupling  effects  as  a  function  of 
subbottom  sediment  type  and  range-rate  of  change  of  sediment  sound  speed. 
The  bulk  of  this  work  is  summarized  in  Ref.  1.  One  outgrowth  of  this 
work  was  that  the  adiabatic  approximation  to  coupled  mode  theory  is 
probably  sufficient  to  address  propagation  problems  concerning  realistic 
lateral  subbottom  variability.  This  clearly  established  the  need  for  the 
capability  to  calculate  propagation  loss  using  the  adiabatic  approximation 
as  a  tool  in  investigating  lateral  variations. 

2.  Sloping  Bottom 

In  the  previous  contract  year  investigations  of  the  importance 

of  mode-mode  coupling  effects  and  how  they  depend  on  bottom  slope  and 

sediment  type  were  begun  and  essentially  completed  (see  Ref.  1).  The 

results  of  this  work  indicated  that  mode  coupling  effects  are  likely  to 

manifest  themselves  for  bottom  slopes  characteristic  of  typical  continental 

rise  and  slope  areas.  This  conclusion  gave  rise  to  an  investigation,  begun 

in  FY78  and  completed  in  FY79,  concerning  the  presence  or  absence  of 

multipath  conversion  effects  in  the  adiabatic  approximation  when  applied 

to  problems  with  sloping  boundaries.  The  intent  of  this  investigation  was 

to  determine  exactly  what  information  was  being  eliminated  by  neglecting 

mode  coupling  effects  in  the  adiabatic  approximation.  Since  this  work 

was  completed  in  FY79  it  is  included  in  this  report.  To  summarize,  at  the 

beginning  of  FY79  the  sloping  bottom  work  had  reached  a  point  where 

(1)  it  was  determined  that  mode  conversion  effects  are  likely  to  be 

important  for  slopes  of  interest,  and  (2)  more  theoretical  study  of 
2  3 

coupled  mode  theory  *  and  its  various  approximations  (adiabatic  and 
boundary  condition  approximations)  with  respect  to  sloping  bottom  problems 
was  called  for. 
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B.  Results  of  FY79  Research 


1.  Lateral  Variability 

In  FY78,  the  need  for  the  capability  to  compute  propagation 
loss  in  the  adiabatic  approximation  was  established;  the  bulk  of  this 
year's  work  has  been  to  build  the  numerical  modeling  tools  necessary  to 
attain  this  capability.  A  result  of  this  work  is  a  numerical  modeling 
tool,  ADIAB,  which  calculates  transmission  loss  in  the  adiabatic  approxi¬ 
mation  for  a  general  waveguide  having  range  variable  bathymetry,  sound 
speed,  attenuation,  etc.  Another  result  has  been  the  selection  of  useful 
power  transport  quantities  for  examining  the  effects  of  lateral  subbottom 
variations.  These  quantities  are  the  power  transported  via  waterborne 
and  bottom  interacting  modes  of  propagation,  and  their  ratios.  Additional 
software  to  be  used  in  conjunction  with  ADIAB  computes  the  power  transport 
quantities  as  a  function  of  range  along  with  additional  modeling  components 
that  permit  estimates  of  the  validity  of  the  adiabatic  approximation.  The 
modeling  capabilities  developed  in  FY79  will  be  used  extensively  in  the 
work  on  lateral  variability  and  sloping  bottoms  to  be  carried  out  in  FY80. 

2 .  Sloping  Bottom 

The  FY79  work  on  sloping  bottom  problems  has  focused  on  a 
theoretical  examination  of  the  mode  coupling  process  as  it  occurs  in 
acoustic  propagation  over  a  sloping  bottom.  The  two  issues  that  were 
examined  are:  (1)  when  one  employs  the  adiabatic  approximation  in  sloping 
bottom  propagation  problems,  is  multipath  conversion  present  or  has  it 
been  eliminated  with  the  neglect  of  mode  coupling?  and  (2)  what  are  the 
implications  of  a  boundary  condition  approximation  inherent  to  coupled 
mode  theory  and  can  the  induced  errors  be  corrected  for? 

Concerning  issue  1,  it  has  been  determined  that  the  adiabatic 
approximation  does  indeed  contain  multipath  conversion  effects.  This 
conclusion  was  reached  through  the  use  and  extension  of  the  work  by  Tindle 
and  Guthrie^*  on  the  analogy  between  ray  theory  and  mode  theory.  This 
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conclusion  is  significant  in  that  the  adiabatic  approximation  to  coupled 
mode  theory  is  substantially  easier  to  implement  than  the  full  mode 
coupling  formulation.  It  is  clearly  important  therefore  to  understand 
the  implications  of  the  approximation  so  that  it  might  be  exploited  when¬ 
ever  valid  and  appropriate.  The  work  concerning  issue  1  is  summarized  in 
Ref.  5. 


Regarding  issue  2,  it  was  determined  that  the  result  of  the 

boundary  condition  approximation  inherent  to  coupled  mode  theory  in 

sloping  bottom  applications  is  that  the  field  given  by  conventional 
2  3 

coupled  mode  theory  ’  does  not  conserve  energy.  The  boundary  condition 
approximation  in  question  involves  the  substitution  of  continuity  of  the 
z  (depth)  derivative  of  the  acoustic  field  for  the  physically  correct 
boundary  condition  of  continuity  of  the  normal  derivative  of  the  field 
at  sloping  interfaces.  This  boundary  condition  approximation  was  thought 
to  be  necessary  to  carry  out  the  partial  separation  of  depth  and  range 
variables  employed  in  the  mathematical  formulation  of  the  theory.  A 
significant  conclusion  of  the  work  reported  here  is  that  if  one  wishes 
to  use  coupled  mode  theory  on  sloping  bottom  problems,  the  adiabatic 
approximation  is  the  only  consistent  form  of  the  theory  that  can  be  used. 
The  adiabatic  approximation  is  consistent  because  the  boundary  condition 
approximation  is  of  the  same  order  as  the  adiabatic  approximation,  whereas 
the  inclusion  of  mode  coupling  effects  is  not.  Therefore,  if  one  wishes 
to  include  mode  coupling  effects,  the  original  formulation  of  coupled 
mode  theory  must  be  altered.  The  FY79  work  on  sloping  bottoms  concluded 
with  the  development  of  a  corrected  coupled  mode  theory  which  includes 
mode  coupling  effects  consistent  with  the  boundary  conditions  satisfied 
by  the  field.  The  FY79  work  on  issue  2  is  summarized  in  Ref.  6  (see  also 
Appendix  A) .  Reference  7  contains  a  more  detailed  account  of  the  work  on 
the  boundary  condition  approximation  as  well  as  more  detailed  versions  of 
the  work  summarized  in  Refs.  1  and  5. 
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c. 


Future  Directions 


With  regard  to  lateral  subbottom  variability,  the  procedure  of  range 
averaging  geoacoustic  parameters  should  be  investigated.  This  investiga¬ 
tion  should  explore  the  relation  between  a  range  averaged  geoacoustic 
bottom  model  and  range  averaging  the  acoustic  field.  Another  important 
problem  concerns  the  relation  between  bottom  loss  obtained  from  multibounce 
measurements  over  a  laterally  variable  ocean  bottom  and  that  obtained  from 
a  range  averaged  geoacoustic  model. 

A  problem  area  involving  sloping  bottoms  is  the  partitioning  of 
energy  for  upslope  and  downslope  propagation  over  continental  slope  type 
bathymetries .  The  change  in  the  amount  of  bottom  interacting  and  water¬ 
borne  energy  as  a  function  of  bottom  slope,  frequency,  and  bottom 
composition  needs  to  be  determined.  In  particular,  one  could  investigate 
(1)  whether  acoustic  energy  traveling  through  the  bottom  and  "popping  out" 
on  the  slope  is  energetically  feasible,  and  (2)  the  relative  amounts  of 
energy  reaching  an  array  position  via  water  and  bottom  interacting  paths 
as  a  function  of  the  "geoacoustics"  of  the  slope. 
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III.  BOTTOM  ROUGHNESS  STUDIES 


The  problem  of  accounting  for  the  effects  of  bottom  roughness  on 
sound  propagation  is  a  very  difficult  one  that  has  received  the  attention 
of  many  researchers  over  a  period  of  years.  Since  the  inception  of  the 
ARL:UT  bottom  interaction  program  a  small  continuous  effort  has  been 
devoted  to  this  specific  problem.  The  goals  of  this  work  have  been  to 
assess  the  usefulness  of  existing  methods  for  accounting  for  bottom 
roughness  and  then  to  carry  out  sensitivity  studies  using  the  most 
promising  of  these  methods.  The  goals  of  the  sensitivity  studies  are 
(1)  to  determine  the  required  level  of  detail  of  description  of  subbottom 
roughness,  and  (2)  to  assess  the  impact  of  bottom  roughness  effects  on 
propagation  problems  relating  to  surveillance. 

A.  Review  of  Previous  Work 


Initial  efforts  were  directed  toward  understanding  the  effects  of 
roughness  at  the  water-sediment  interface  within  the  context  of  ARL:UT 
bottom  loss  and  normal  mode  models,  which  treat  the  sediment  as  a  fluid. 

The  best  available  method  for  treating  a  randomly  rough  two-fluid  inter- 

13 

face  was  found  to  be  a  smoothed  or  averaged  boundary  condition  approach. 

3 

This  approach  was  implemented  in  FY77;  sensitivity  studies  suggested  that 
the  roughness  of  the  water-seaiment  interface  is  not  important  at  low 
frequencies  for  high  porosity  sediments  (clays  and  silts)  such  as  those 
found  in  deep  ocean  basins  (abyssal  plains) . 

In  FY78  it  was  suggested  that,  for  regions  of  thin  sediment  cover, 
such  as  the  Northeast  Pacific,  the  most  important  roughness  effects  might 
occur  at  the  sediment-substrate  interface.  Application  of  the  two-fluid 
theory  was  not  attempted  because  it  was  recognized  that  the  solid  (shear 
wave)  properties  of  the  substrate  were  important.  Work  to  develop  an 
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average  boundary  condition  approach  to  describe  a  randomly  rough  fluid- 
solid  interface  was  initiated  and  completed  to  the  stage  of  equations 

4 

suitable  for  numerical  computation. 

B.  Results  of  FY79  Research 


Computer  software  to  implement  the  average  boundary  condition 
approach  to  a  randomly  rough  fluid-solid  interface  was  developed  and 
tested.  The  computer  program  calculates  the  reflection  and  transmission 
coefficients  of  the  coherent  component  of  the  sound  field.  Techniques 
were  also  developed  to  use  these  coefficients  to  modify  the  ARL:UT  fluid 
sediment  bottom  reflection  loss  model"*  to  include  the  effects  of  substrate 
roughness.  Results  of  the  sensitivity  studies  carried  out  under  this  task 
will  appear  in  a  future  publication. 
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IV.  BOTTOM  INTERACTION  SENSITIVITY  STUDIES 


A  continuing  aspect  of  the  bottom  interaction  studies  at  ARL:UT 
has  been  the  use  of  various  "measures"  of  bottom  interaction  as  vehicles 
to  test  the  importance  of  one  or  more  parameter  variations  or  uncertainties. 
Initially  a  commonly  used  measure  was  propagation  loss,  and  work  was 
carried  out  to  determine  the  sensitivity  of  propagation  loss  to  variations 
in  bottom  loss.  More  recent  studies  have  explored  the  acoustical  impor¬ 
tance  of  various  sea  floor  parameters  such  as  density,  sound  speed,  and 
absorption  gradient  using  bottom  reflection  loss  as  a  measure  of  bottom 
interaction . 

In  FY76  a  study  of  the  sensitivity  of  bottom  reflection  loss  to 
subbottom  parameter  variations  began.  Initially  the  work  concentrated  on 
the  properties  of  a  fluid  sediment  and  was  carried  out  using  a  sophisti¬ 
cated  bottom  reflection  loss  model1  developed  at  ARL:UT  for  this  purpose, 

2  3 

Various  sea  floor  parameters  such  as  density  gradient,  sound  speed  and 

4  3  6 

absorption  gradients,  and  substrate  rigidity  *  have  been  successfully 
studied . 

A.  Review  of  Previous  Work 


In  FY78  the  direction  of  these  studies  turned  toward  the  inclusion  of 
sediment  rigidity  (shear  waves)  in  the  bottom  reflection  loss  model.  The 
necessary  theoretical  work  to  develop  a  numerical  bottom  reflection  loss 
model  for  use  in  studying  the  acoustic  properties  of  solid  sediments  was 
initiated,  completed,  and  implemented . ^  The  assumptions  underlying  this 
model  restrict  its  validity  to  frequencies  above  about  10  Hz.  Theoretical 
work  to  extend  the  frequency  range  to  lower  frequencies  was  also  completed. 
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IniCial  studies  of  a  hypothetical  turbidite  layer  showed  that  sediment 
shear  waves  were  more  important  than  heretofore  thought,  particularly  for 
regions  of  thin  sediment  cover.  Preliminary  results  suggested  that  the 
major  mechanism  for  sediment  shear  wave  excitation  is  compressional  wave 
conversion  at  the  sediment-substrate  interface  and  that  the  water -sediment 
interface  can  be  accurately  represented  as  a  fluid-fluid  interface. 

B .  Results  of  FY79  Research 

During  FY79  computer  software  was  developed  to  be  used  in  preliminary 
studies  of  bottom  interaction  effects  on  coherence.  The  procedure  followed 
uses  the  ARL:UT  normal  mode  model  to  calculate  the  sound  field  for  realistic 
environmental  and  subbottom  acoustic  properties.  The  coherence  between  two 
receivers  at  the  same  range  but  different  depths  is  then  computed.  The 
averaging  procedure  simulates  a  time  average  of  the  received  signal. 

Doppler  effects  are  not  included  at  this  time  so  that  multipath  and  bottom 
interaction  effects  could  be  isolated. 


The  major  effort  in  FY79  concentrated  on  the  effect  of  sediment 
rigidity  on  bottom  reflection  loss.  This  work  is  reported  in  detail  in 
Ref.  8,  which  is  being  segmented  into  several  smaller  reports  (Refs.  9-12) 
for  journal  publication.  References  9  and  10  appear  as  Appendices  B  and  C 
of  this  report.  The  results  of  FY79  work  are  summarized  below. 

1 .  The  Effect  of  Sediment  Rigidity  on  Bottom  Reflection  Loss  in 
a  Typical  Deep  Sea  Sediment 

A  study  of  the  effect  of  sediment  rigidity  on  bottom  reflection 
loss  (RL)  from  typical  deep  sea  sediment  types  shows  that  sediment  shear  (S) 
wave  excitation  is  important  for  thin  sediment  layers  but  is  negligible  for 
thick  layers.  The  major  mechanism  for  S  wave  excitation  is  compressional 
(P)  wave  conversion  at  the  sediment-substrate  interface.  Little  energy  is 
coupled  into  S  waves  at  the  water-sediment  interface. 


The  mechanism  for  S  wave  excitation  provides  a  means  for 
quantifying  the  categories  of  thick  and  thin  layers.  It  is  the  amplitude 
of  the  P  wave  at  the  substrate  interface  that  matters.  If  the  P  wave 
amplitude  at  the  substrate  interface  is  significant,  the  layer  is  thin; 
if  it  is  small,  the  layer  is  thick.  The  categories  of  thick  and  thin,  then, 
will  depend  upon  the  wave  frequency  (through  P  wave  attenuation),  grazing 
angle  (through  the  depth  of  the  turning  point),  and  the  sediment  thickness. 

For  physically  thin  sediment  layers  the  dependence  of  RL  on 
frequency  (at  a  given  grazing  angle)  exhibits  high  and  low  frequency 
regimes  with  strikingly  different  properties.  The  existence  of  these 
frequency  regimes  is  the  result  of  interference  effects  related  to  S  wave 
propagation  within  the  sediment.  At  low  frequencies  sediment  shear  wave 
attenuation  is  small  and  interference  effects  dominate  the  dependence  of 
RL  on  frequency.  Large  peaks  in  RL  (=20  dB)  occur  with  a  separation  in 
frequency  related  to  changes  in  S  wave  phase  of  2ir.  At  higher  frequencies 
the  increased  shear  wave  attenuation  decreases  the  magnitude  of  this  peak 
structure.  At  high  frequencies  the  oscillatory  structure  disappears  as 
S  waves  excited  at  the  sediment-substrate  interface  are  totally  absorbed 
within  the  sediment,  resulting  in  an  almost  constant  additional  loss 
compared  to  the  RL  of  a  fluid  sediment. 

The  total  attenuation  of  the  S  wave  traveling  one  way  through 
the  sediment  provides  a  quantitative  means  of  separating  these  frequency 
regimes.  For  total  S  wave  attenuation  greater  than  20  dB  an  S  wave 
generated  at  the  substrate  interface  is  essentially  absorbed  within  the 
sediment  before  striking  the  other  interface.  Since  attenuation  is 
proportional  to  frequency,  this  defines"’  a  high  frequency  regime  in  which 
the  propagation  of  sediment  S  waves  can  be  neglected;  however,  their 
excitation  is  still  an  important  loss  mechanism.  At  low  frequencies  S 
wave  attenuation  is  less  than  20  dB  across  the  layer  and  the  additional 
interference  effects  due  to  the  propagation  of  S  waves  through  the 
sediment  are  important. 


2. 


A  Ray  Path  Analysis  of  Sediment  Shear  Wave  Effects  on  Bottom 
Reflection  Loss 


A  theoretical  treatment  of  the  reflection  and  transmission 
coefficients  at  the  water  and  substrate  interfaces  by  means  of  an  expansion 
in  the  small  sediment  S  wave  velocity  provides  a  basis  for  developing  a 
ray  path  model  of  the  effect  of  sediment  rigidity  on  bottom  reflection 
loss.  The  basic  mechanism  for  S  wave  excitation  emerges  as  P  wave 
conversion  at  the  substrate  interface.  The  dependence  on  S  wave  velocity 
and  attenuation,  and  P  wave  attenuation,  are  predicted  by  the  ray  model 
and  verified  in  detail  by  results  obtained  from  the  computational  bottom 
reflection  loss  model.7 

The  bottom  reflection  loss  is  found  to  be  due  to  the  interference 

of  three  waves  in  the  water:  the  Aq  wave,  similar  to  that  generated  by  a 

fluid  sediment,  the  A^  wave,  whose  amplitude  depends  on  the  S  wave  velocity 

to  the  first  power,  and  the  smaller  A^  wave,  whose  amplitude  depends  on  the 

S  wave  velocity  to  the  second  power.  The  interference  between  the  Aq  and 

A^  waves  produces  the  large  interference  peaks  in  the  low  frequency  regime 

of  RL  noted  above.  The  A^  wave  provides  an  additional  modulation  to  the 

peak  structure.  At  high  frequencies  the  A^  and  A^  waves  become  negligible 

as  the  attenuation  of  the  shear  wave  increases.  The  A  wave,  however, 

o 

remains  and  carries  the  effect  of  shear  waves  excited  at  the  substrate 

and  absorbed  within  the  sediment.  The  A  wave  describes  the  increased 

o 

loss  in  the  high  frequency  regime  of  RL. 

3 .  Sensitivity  Studies 


Studies  to  identify  significant  parameters  determining  bottom 
reflection  loss  in  thin  sediment  layers  were  completed  in  FY79.  For 
homogeneous  sediment  layers,  at  high  frequencies,  the  sediment  S  wave  and 
P  wave  velocities,  sediment  P  wave  attenuation,  and  substrate  P  and  S  wave 
velocities  are  important  parameters.  Varying  these  parameters  through 
their  geophysically  allowed  ranges  produces  a  change  in  bottom  reflection 
loss  of  more  than  0.5  dB  relative  to  that  calculated  using  the  value  at 


the  center  of  their  range  attenuation.  At  lower  frequencies  S  wave 
propagation  within  the  sediment  becomes  important  and  the  sediment  S 
wave  attenuation  also  becomes  an  important  parameter.  Studies  of  thin 
sediment  layers  with  gradients  show  that  gradients  of  compress ional  and 
shear  wave  velocity  and  shear  wave  attenuation  are  also  important  param¬ 
eters  . 

Studies  indicating  the  accuracy  required  in  specifying  the  geoacoustic 
parameters  of  thin  sediment  layers  were  also  completed.  Bottom  reflection 
loss  from  solid  sediments  was  found  to  be  particularly  sensitive  to  water 
sound  speed,  sediment  thickness,  sediment  surficial  density,  sediment 
surficial  compressional  wave  speed,  and  substrate  shear  velocity.  Frac¬ 
tional  changes  of  these  parameters  induced  much  larger  fractional  changes 
in  the  magnitude  of  the  reflection  coefficient  of  the  ocean  bottom. 

Several  parameters  induced  negligible  changes  in  bottom  reflection  loss 
when  varied  over  their  geophysically  expected  ranges.  Thes“  very  low 
sensitivity  parameters  were:  all  attenuations,  the  gradients  of  all 
attenuations,  the  gradients  of  sediment  density,  and  the  gradient  of  the 
sediment  compressional  wave  velocity.  The  remaining  parameters  fell 
between  these  limits  of  sensitivity,  producing  fractional  changes  in  the 
reflection  coefficient  of  the  same  size  as  the  fractional  change  in  the 
parameter. 

4 .  Effect  of  Sediment  Rigidity  on  Interface  Wave  Excitation 

Sediment  rigidity  was  also  found  to  influence  the  energy  lost  to 
interface  waves.  The  Stoneley  wave  peak,  generated  in  the  fluid  sediment 
case,  merges  into  the  generally  increased  bottom  reflection  loss  in  the 
solid  sediment  case.  The  characteristic  peak  in  the  dependence  of  bottom 
reflection  loss  on  grazing  angle,  seen  in  the  fluid  sediment  case,  cannot 
be  unambiguously  identified  in  the  case  of  a  solid  sediment  with  realistic 
geophysical  parameter  values.  This  work  goes  some  way  toward  explaining 
why  the  previously  predicted  Stoneley  wave  reflection  loss  peaks  have  not 
been  observed,  and  to  clarify  the  true  role  of  these  interface  waves  in 
the  reflection  loss  process. 
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A  consistent  coupled  mode  theory  of  uound  propagation  for  a  class  of  non- 
separable  problems*^ 

Steven  R.  Rutherford  and  Kenneth  Z.  Hawker 
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This  article  examines  the  effects  of  boundary  condition  approximations 
that  arise  whenever  the  coupled  mode  theories  of  fierce  and  Milder  are 
applied  to  propagation  problems  involving  range  dependent  boundaries. 

This  boundary  condition  approximation  requires  that  the  depth  derivative, 
rather  than  the  normal  derivative  of  the  field,  be  continuous  across  a 
sloping  interface.  The  approximation  is  necessary  in  order  to  carry  out 
the  partial  separation  of  depth  and  range  variables  effected  in  the  mathematical 
formulation  of  the  theory.  This  article  will  show  that  a  consequence  of 
this  approximation  is  that  conventional  coupled  mode  theory  applied  to 
dissipationless  media  with  nonhorizontal  boundaries  does  not  conserve  energy. 

It  is  shown  that  a  correction  to  coupled  mode  theory  can  be  derived  such 
that  the  proper  boundary  conditions  are  satisfied  and  energy  is  conserved 
to  first  order  in  the  local  bottom  slope.  Moreover,  the  corrections  are 
not  prohibitive  in  terms  of  added  computational  complexity.  Numerical 
examples  are  presented  which  illustrate  the  nonconservation  of  energy  effect 
and  the  corrections  to  the  theory. 

PACS  numbers : 
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INTRODUCTION 


In  the  coupled  mode  theory  of  underwater  sound  propagation  proposed 
1  .  2 

independently  by  Fierce  and  Milder,  a  certain  form  of  the  boundary  conditions 
satisfied  by  the  field  was  assumed,  namely,  continuity  of  pressure  and 
vertical  component  of  particle  velocity.  This  assumption  allowed  a  partial 
separation  of  the  depth  and  radial  variables  and  resulted  in  a  workable 
theory.  In  range  dependent  media  having  horizontal  boundaries  and  interfaces, 
the  assumed  form  of  the  boundary  conditions  was  physically  correct;  however, 
when  applied  to  problems  having  range  variable  boundaries  and  interfaces, 
the  assumed  boundary  conditions  were  approximations  to  the  proper  ones 
to  be  satisfied  by  the  field  at  fluid-fluid  interfaces  (continuity  of  normal 
stress  and  displacement).  The  intent  of  this  article  is  to  explore  the 
consequences  of  this  approximation  to  the  physically  correct  boundary  conditions. 
It  is  shown  that  one  consequence  of  the  approximation  is  nonconservation 
of  energy  within  the  original  formalism  of  coupled  mode  theory.  It  is 
also  shown  that  a  consistent  correction  to  the  theory  can  be  made  such 
that  the  field  so  obtained  satisfies  both  the  proper  boundary  conditions 
and  conserves  energy  to  first  order  in  the  slopes  of  nonhorizontal  boundaries 
and  interfaces  that  are  present  in  the  problem. 

3-8 

Of  the  recent  applications  of  coupled  mode  theory,  most  involved 
media  with  horizontal  boundaries;  hence,  as  shall  be  shown,  the  issue  of 
nonconservation  of  energy  did  not  arise.  The  adiabatic  approximation  to 
coupled  mode  theory  (see  Refs.  9-12)  has  been  applied  to  many  types  of 
acoustic  propagation  problems.  The  nonconservation  of  energy  effect  does 
not  arise  in  adiabatic  mode  theory  (see  Ref.  9)  because  the  assumed,  approximate 
boundary  conditions  are  consistent  with  the  assumption  of  adiabacity. 


I. 


BACKGROUND 


In  the  coupled  mode  theories  of  Pierce  and  Milder,  the  acoustic  velocity 
potential  in  a  range  dependent,  fluid  medium  is  expressed  as 


'J'(r.z)  a£Rn(r)  4>n(z;r) 

n 

with  Rn(r)  and  <$>n(z;r)  satisfying 
2 
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*n(z;r)  =  o.  (3) 

The  summation  in  Eq.  (1)  extends  over  all  discrete  modes.  In  Eqs.  (1) 
through  (3)  a  cylindrical  coordinate  system  has  been  assumed  and  the  azimuthal 
symmetry  of  the  medium  used  to  reduce  the  problem  to  two  dimensions.  The 
velocity  potential  ^  satisfies  the  following  partial  differential  equation 


V2<J>  ♦  k2(r,z)  =  0  (4) 

with  k(r,z)  being  the  wave  number  of  the  medium  (k(r,z)»<*Vc(r,z)) .  As 

indicated  by  Eq.  (4)  the  developments  of  this  article  will  consider  acoustic 

propagation  in  a  source  free  region  of  space.  The  extension  to  include 

a  source  is  straightforward.  The  boundary  conditions  on  <t>n  are  that  P<J>n 

and  ^  /3zbe  continuous  across  any  interfaces.  The  functions  <J>  are  assumed 
n  3  n 

to  form  an  orthonormal,  complete  set  of  functions  such  that 
co 


/N>  <J>  (z; r)  dz  =  6 

n  m  n,m 
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o 

where  P  is  the  density  of  the  medium,  assumed  to  have  step  function  dis¬ 
continuities  at  boundary  interfaces.  If  the  eigenfunctions  corresponding 
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Co  Che  discreCe  spectrum  do  noc  fora  a  complete  set  (due  Co  Che  presence 

of  a  continuous  spectrum),  Che  results  throughout  this  paper  remain  correct 

but  only  approximate.  The  A  and  B  in  Eq.  (2)  are  range-dependent  coupling 

mn  inn 

coefficient  matrices  defined  by 


B  (r) 
mnv  J 


A  (r)  = 
mnk  J 


n  . 

!Tdz 

(6) 

00 

„  2 

P<J> 

in 

n  , 

— z-  dz  . 

(7) 

Implicit  in  the  partial  separation  of  variables  described  by  Eqs.  (1) 
through  (4)  is  that  the  field  ^  satisfies  boundary  conditions  such  that 
P>P  and  Sty  are  continuous  throughout  the  waveguide.  These  boundary  conditions 
are  equivalent  to  requiring  that  the  pressure  and  z-coraponent  of  particle 
velocity  be  continuous  functions  of  depth  and  range.  These  boundary  conditions 
are  the  physically  correct  ones  for  range  dependent  waveguides  having  horizontally 
stratified  boundaries.  However,  if  one  wishes  to  consider  nonstratif ied 
boundaries,  the  correct  conditions  for  ^  to  satisfy  at  surfaces  of  discontinuity 
are  continuity  of  Pty  and  8^/8n,  where  3/3n  is  the  normal  derivative  operator 
given  by 


JL 

3n 


J i  +  A 


=  —  -  H(r)  — 


H(r) 


l9z 


3r 


(8) 


In  Eq.  (8),  H(r)  is  the  local  slope  of  a  surface  of  discontinuity  defined 
by  z**H(r)  with  H(r)“dH/dr;  H(r)  might,  for  example,  be  the  water  depth 
at  range  r. 

Of  course,  if  is  required  to  satisfy  boundary  conditions  given  by 
Eq.  (8),  the  partial  separation  of  the  range  and  depth  variables  implied 
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by  Eq.  (1)  is  not  possible.  In  this  situation  the  variables  are  not  separable 
and  an  exact  solution  to  Eq.  (4)  having  the  assumed  form  of  'I'  given  in 
Eq.  (1)  cannot  be  obtained.  A  natural  approximation  to  make  when  applying 


coupled  mode  theory  to  problems  involving  range  dependent  boundaries  is 
to  approximate  the  normal  derivative  of  the  field  at  the  boundary  with 


by  Eqs.  (1)  through  (4)  can  be  applied  to  problems  having  nonhorizontal 


boundaries.  It  is  usually  reasoned  that,  for  small  boundary  slopes,  the 
errors  incurred  from  the  approximation  of  Eq.  (9)  will  be  negligible. 

It  is  the  purpose  of  this  paper  to  explore  the  implications  of  the  approximation 
of  Eq.  (9)  and  to  show  that  a  correction  to  coupled  mode  theory  must  be 


applied  whenever  nonhorizontal  boundaries  are  present. 

The  coupling  coefficients  of  Eqs.  (6)  and  (7)  are  quantities  that 

depend  on  the  radial  rates  of  change  of  the  medium.  It  will  be  shown  that 

•  *2 

B  and  A  have  components  of  order  H  and  H  ,  respectively.  Therefore, 
mn  ran 

the  inconsistency  of  the  approximation  in  Eq.  (9)  is  immediately  apparent 
since  terms  of  order  H  and  higher  are  being  neglected  in  Eq.  (9)  while 
they  are  being  retained  in  Eq.  (2)  with  the  coupling  coefficients.  An 
implication  of  this  inconsistency  is  that  if  one  insists  on  satisfying 
the  boundary  conditions  on  'I*  only  to  zero  order  in  H,  as  given  by  Eq.  (9), 
then  the  only  consistent  treatment  of  the  radial  equation  is  one  in  which 
mode  coupling  effects  are  ignored  (i.e.,  the  adiabatic  approximation). 

In  the  next  section  of  this  article  it  will  be  shown  that  a  consequence 


of  the  inconsistent  use  of  Eq.  (9)  is  a  nonconservation  of  energy  within 
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the  conventional  coupled  mode  theory  formalism.  This  inconsistency  problem 
has  not  arisen  until  now  because  previous  applications  of  coupled  mode 
theory  were  ones  in  which  boundary  conditions  involving  the  normal  derivative 
either  did  not  appear,  or  else  the  subsequent  assumption  of  adiabatic  propagation 
avoided  the  question. 

II.  CONSIDERATIONS  OF  ENERGY  TRANSPORT 

In  this  section  we  shall  examine  energy  transport,  and  the  question 
of  energy  conservation  in  dissipationless  media,  using  the  form  of  coupled 
mode  theory  outlined  in  the  previous  section. 

It  is  of  course  recognized  that  preservation  of  an  intact  fluid-fluid 
interface  requires  continuity  of  both  displacement  and  normal  stress,  leading, 
in  the  time  harmonic  problem  considered  here,  to  continuity  of  d^/dn  and 
'I'.  In  this  sense,  then,  it  is  required  that  dip/d n  be  continuous  and  it 
should  not  be  surprising  that  inconsistencies  develop  when  this  requirement 
is  not  met.  However,  in  as  much  as  coupled  mode  theory  has  been  developed 
for  the  approximate  continuity  condition,  it  becomes  necessary  to  explore 
the  consequences  of  this  requirement.  Moreover,  it  is  not  a  priori  obvious 
to  what  order  in  H(r),  the  local  interface  slope,  inconsistencies  may  develop 
or  whether  the  inclusion  of  mode  coupling  at  all  is  consistent  with  the 
approximation  ^/^n  =  3'i'/3z. 

Following  Ref.  9,  we  consider  the  time  averaged  (rms)  radial  energy 
flux  Jf  given  by^ 

Jr  =  j  Re(pv*)  =  -  j  wp  .  (10) 

where  v^  is  the  radial  velocity  and  p  the  pressure.  Also  of  interest  will 
be  the  total  power  transported  through  a  cylindrical  surface  of  radius 


A3 


r,  centered  at  the  origin  and  extending  from  the  sea  surface  (z=0)  to  infinity. 
This  quantity  P  is  given  by 


2ir  oo 


P  =/  d df  dz  rj  =  2*rf  dz  J  (r,z) 

Oo  -u  a 


(ii) 


where  it  must  be  remembered  that  integrals  of  the  form  /  dz(*)  are  abbreviated 

O  it  .. 


no 


tation  for  a  sura  of  integrals  over  the  contiguous  layers  fQ  dz( • )  +  /^dz(  • ) , 
Upon  differentiating  with  respect  to  r  and  using  Eq.  (10)  we  have  for  the 


single  interface  (water-sediment)  problem 


3P 

57 


=  -ttui  imj/  dz  p  ^ 


3^* 

3r 


+  r 


3fl> 


3r 


+ 


34 ' 

3r2 


li 


(12) 


_  nil 

where  H  and  H  refer  to  limits  approaching  the  water-sediment  interface 

(z=H(r))  from  above  and  below,  respectively.  The  wave  equation  can  now 

2  2 

be  used  to  substitute  for  3  i//*/3r  ,  and  after  routine  manipulations  one 


obtains 


3P 

57 


=  -irwrH  -  'i'*)] 


(13) 


+TTU)r  ImJ  J  dz  ^p^— , 

with  H=dH/dr  and  =  (3tfy3r)|  ,  etc. 

In  passing  from  Eq.  (12)  to  Eq.  (13)  it  has  been  assumed  that  is 
continuous  across  the  interface  z=H(r).  Without  further  assumptions,  for 
the  moment,  as  to  the  continuity  of  d'P/dz,  one  has,  from  Eq.  (13)  upon 
integration  by  parts  in  the  second  term,  the  result 

=  -TTwrlmlp^jHO^  -  (J/*)  -  0J»L'  -  >Jy)j| 

where  ^  =  [34<z ,  r ) /  ^z  ]  ,  etc . 

z+H 

If,  as  has  been  assumed  in  conventional  coupled  mode  theory,  3i///3z 
is  taken  to  be  continuous  at  all  interfaces,  the  last  term  in  Eq.  (14)  vanishes 


(14) 
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and  one  has  the  result, 

BP  I 

37=  -™r  ImJpj^j  Hflpj  -  i|>2)|  .  (15) 

At  this  stage  it  will  be  useful  to  introduce  the  partially  separated  solution 
given  by  Eq.  (1).  After  carrying  out  some  straightforward  manipulations, 
as  described  in  Appendix  A,  one  obtains 


SP 

3r 


nm , 


r  Q  (r) 


+  r  B  P  (r) 
nm  nm 


(16) 


where  the  superscripts  A  and  S  refer  to  the  antisymmetric  and  symmetric 

parts  of  the  matrices.  The  range  dependent  factors  Q  and  P  are  defined 

nm  nm 

V  ■  2“»IW  (17) 


as 


P 

nm 


=  2iIm(R  R* 
n  m 


+  R  R  ) 
m  n 


(18) 


while,  as  shown  in  Appendix  A,  the  matrices  in  Eq.  (16)  are  given  by 


(B  -  A  )A  =  -  Up 
nm  nm'  2 


n  m  m  '  m  n 


;(2). 


BS  .  >  Hp,(l  -  ♦ 

nm  2  pl\  pj  n  Y 


O)  *0) 

m 


(19) 


(20) 


for  the  waveguide  geometry  of  Fig.  1. 

In  a  general  problem  where  H^O  there  is  clearly  no  reason  to  expect 
the  right-hand  side  of  Eq.  (16)  to  vanish.  Thus,  in  general,  BP/Sr^O  and 
the  (radial)  energy  flux  passing  through  a  semiinfinite  cylindrical  surface 
is  not  independent  of  range.  Of  course,  in  the  horizontally  stratified 
case  H  vanishes  identically  and  energy  is  properly  conserved.  Similarly, 
in  the  adiabatic  approximation  to  a  nonhorizontally  stratified  problem, 
all  mode  coupling  terms  are  ignored  and  once  again  energy  is  conserved. 
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Even  Chough  we  have  not  yet  presented  a  coupled  mode  theory  which 
satisfies  the  proper  continuity  relations  at  a  fluid-fluid  interface,  the 
presence  of  (only)  boundary  terms  on  the  right-hand  side  of  Eq.  (16)  leads 
one  to  inquire  as  to  the  respective  roles  of  the  continuity  of  3^/9z  and  9^/3 n 
at  the  boundary.  Returning  to  Eq.  (14)  and  using  Eq.  (8)  one  has. 


hoK  -  i1*)  -  (r*  -  «;*)  =  2 

i  t  i  Z  on 


'*  Ih* 


5T-h- 


Thus,  when  9*JV  is  assumed  continuous  at  z=H(r),  we  have  at  once  that 
3P/&r=0  and  energy  is  conserved.  The  motivation  for  developing  a  coupled 
mode  theory  based  on  continuity  of  3^/3n  rather  than  3 3z  is  therefore 
clearly  established. 

As  a  final  point  we  shall  clarify  the  role  of  the  water-sediment  interface, 
and  the  continuity  of  3^/Sz  imposed  there,  in  determining  the  overall  flow 
of  acoustic  energy.  It  is  easily  verified  by  direct  computation  that 
where  is  the  time  averaged  energy  flux  vector,  the  radial  component  of 
which  is  given  by  Eq .  (10).  Hence,  by  Gauss'  theorem  the  surface  integral 
of  the  normal  component  of  ,£,  integrated  over  a  closed  surface,  vanishes. 

We  now  choose  as  the  surface  a  cylinder  of  radius  r,  centered  at  the  origin, 
and  extending  from  the  surface  (z=0)  to  infinity.  Since  only  at  the  boundary 
z=H(r)  can  any  component  of  J  be  discontinuous,  and  because  J  vanishes 
at  z=0,  and  as  z"*”,  we  have 


n 

/  d6  J  dz  rJ  +  f|jl  -  Ji  1  dA  =  0 

n  a  J  C  *  * 


o  o 


where  Jx  is  the  energy  flux  normal  to  the  surface  defined  by  the  water- 
sediment  interface,  and  the  ^ ^ superscripts  refer  to  the  fluxes  normal 
to  the  interface  as  approached  from  above  and  below. 


_  <  i mm*  , 


After  straightforward  computation,  the  net  energy  flux  across  the  water- 
sediment  interface,  j^-j£  ,  is  found  to  be 

Jj.  -  J*  =  -  A  \A  0  ♦  BS  P  I  (271 

~  4/uirfcJ'nm  nm'  ™  nm  'nm  |  •  U?J 

Upon  introducing  Eq.  (27)  into  Eq.  (24),  the  expression  for  the  radial 
rate  of  change  of  total  power  transmitted  through  the  cylindrical  surface 
must  obey  the  conservation  equation, 

9P  iCJTT  t— >1/  '  \A  Si 

57  '  nB  '  A"")  r  Q™  *  rB™  Pnm  I  '  °  '  (28) 

However,  as  previously  shown  in  Eq.  (16),  the  second  term  in  Eq.  (28)  is 
simply  - 3p/ 3r  so  the  total  energy  balance  equation  is  satisfied  identically. 

It  is  appreciated  that  the  conventional  coupled  mode  theory  outlined 
in  Section  I,  the  energy  considerations  given  here,  and  the  consistent 
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coupled  mode  theory  to  be  developed  in  Section  II,  have  all  been  restricted 
to  lossless  media.  In  practice  when  significant  acoustic  energy  interacts 
with  the  sea  floor,  the  volume  absorption  within  the  sediments  must  be 
taken  into  account.  In  horizontally  stratified  situations  this  is  accomplished 
most  easily  by  letting  wave  numbers  take  on  an  imaginary  part,  k"*k+ia, 
and  obtaining  the  corresponding  imaginary  correction  to  the  eigenvalues, 
k^+ i6n,  by  first  order  perturbation  theory.  It  is  not  difficult  to 
convince  oneself  that  a  similar  result  will  hold  for  a  coupled  mode  theory 
when  absorption  is  small.  This  subject  will  be  addressed  in  a  later  paper. 

In  summary,  we  have  shown  in  this  section  that  (1)  conventional  coupled 
mode  theory,  employing  continuity  of  <W/3z  at  all  interfaces,  does  not 
properly  conserve  total  radial  energy  flux  in  a  dissipationless  problem, 

(2)  the  anomalous  energy  loss  (right-hand  side  of  Eq.  (16))  can  be  expressed 
in  terms  of  mode  coupling  coefficients  and  radial  functions  Rn>  together 
with  the  interface  slope  H,  (3)  total  energy  is  conserved  if  one  includes 
energy  lost  into  the  bottom  as  a  result  of  the  application  of  the  incorrect 
boundary  condition,  and  (4)  if  3ify8n  is  assumed  continuous  at  the  interface, 
the  total  radial  energy  flux  is  correctly  conserved. 

Finally,  we  wish  to  remark  that  as  a  consequence  of  the  use  of  continuity 
conditions  on  3^/3z,  resulting  in  a  theory  which  fails  to  conserve  energy 
to  order  H,  when  one  is  concerned  with  small  but  finite  slopes,  it  will 
be  manifestly  inconsistent  to  use  conventional  coupled  mode  theory  (wherein 
the  coupling  coefficients  are  corrections  of  order  H)  except  in  the  adiabatic 
approximation.  In  this  theory,  we  can  already  see  that  the  first  order 
corrections  to  the  adiabatic  theory  will  be  proportional  to  B  and  thence 
at  least  to  H.  Thus  retaining  mode  coupling  as  the  conventional  theory 
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will  entail  retaining  some  terms  of  order  H  while  ignoring  other  terms 
of  the  same  order  (arising  through  the  continuity  conditions  imposed  at 
sloping  interfaces).  Of  course  in  range  variable  problems  not  involving 
sloping  interfaces  this  inconsistency  does  not  arise. 

III.  FIRST  ORDER  CORRECTIONS  TO  COUPLED  MODE  THEORY 

In  this  section  a  first  order  correction  to  coupled  mode  theory  is 
proposed  and  derived.  It  will  be  shown  that  with  an  improvement  on  the 
theories  of  Pierce  and  Milder,  the  field  can  be  obtained  in  terms  of  coupled 
modes  so  that  the  proper  boundary  conditions  are  satisfied  and  energy  is 
conserved  (both  to  first  order  in  the  boundary  slope  H)  and  a  consistent 
theory  is  obtained. 

To  derive  the  appropriate  corrections  to  the  theory,  Eq.  (1)  is  replaced 

Wr,z)  eZ)Xn(r)  $n  (z;r)  .  (29) 

n 

When  V  is  required  to  satisfy  the  proper  boundary  conditions  (continuity 
of  normal  stress  and  displacement)  to  first  order  in  H,  the  appropriate 
conditions  on  and  its  z  derivative  at  a  sloping  interface  defined  by 
z=H(r)  are  found  to  be  that 


pi> 


and 


ft* 


(30) 


n  Xn  n 


be  continuous  across  z=H(r).  In  these  expressions  the  dot  and  prime  symbols 

denote  differentiation  with  respect  to  range  and  depth,  respectively. 

The  second  of  Eqs.  (30)  follows  from  the  fact  that  $  is  first  order  in 

n 

H. 

Now,  assume  that  the  depth  equation  ♦  satisfies  the  following  differential 


equation. 


[i 


*2  +  k2  (r, z)  -  o*  (r) 
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<&n  (z;r)  =  0 


(31) 


m 


Equation  (31)  is  of  the  same  form  as  the  differential  equation  assumed 
in  the  developments  of  Ref.  1,  the  only  difference  being  that  it  satisfies 
the  modified  boundary  conditions  of  Eqs.  (30). 

If  Eqs.  (29)  and  (31)  are  used  in  Eq.  (4),  one  obtains  the  following 


equation 


oh) 

n  n/ 


0 

n 


(32) 


It  is  now  convenient  to  express  0  as 


0  =0  +60 


(33) 


where  0  satisfies  Eq.  (3)  with  0‘  and  p0  continuous  at  z=H(r).  The  60 
n  n  n  n 

is  a  function  of  order  H  that  vanishes  in  the  limit  of  horizontally  stratified 
boundaries.  Equations  (3)  and  (33)  in  Eq.  (31)  yield  a  differential  equation 


for  6<P  . 
n 


(p"2  ■ 

The  boundary  conditions  on  <$0^  are  obtained  by  using  Eq.  (33)  in  Eqs.  (30) 
and  retaining  only  terms  of  order  H.  For  a  nonhorizontal  boundary  defined 
by  z=H(r)  separating  two  fluids  of  densities  and  ,  as  shown  in  Fig.  1, 
the  boundary  conditions  on  60^  are  found  to  be 


(34) 


p,60(1)  =  p960(2)  u.  , 

1  n  K2  rn  z=H(r) 


60 


.(D-  60'<2> 


JLh  p^Vi  - 

X  ^n  v  p-/z=H(r) 


(35) 


n  -  '  -2' 

In  Eq.  (35)  the  superscripts  1  and  2  mean  values  of  the  function  approached 
from  above  and  below  the  interface  at  z=H(r),  respectively. 

Equations  (34)  and  (3)  can  be  combined  to  give 


0  60"  -  60  0"  =  (o2  -  k2)  0  60  +  fo2  -  k2^  0  0 
Ym  \  n  m  /  rm  Yn  \  n  n  /  rn  rm 


(36) 
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If  Eq.  (36)  is  multiplied  by  p  and  integrated  over  depth,  the  resulting 
equation  is 


E 

i 

where 


p.<i>(i)|V(i)  -  6<r(i+1)] 

KiYm  [  Yn  n  j 


=  (a2  -  k2)  I  +  (o2  -  k2)6 
n  m  mn  n  n  n,m 


z=H. 
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mn 
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pii  6<b  dz 
Yn 


(37) 


(38) 


In  Eq.  (37)  the  summation  is  over  all  sloping  interfaces  in  the  waveguide 
defined  by  z“H£(r).  The  superscripts  i  and  i+1  denote  values  of  the  functions 
approached  from  above  and  below  the  interface  z=H^(r).  If  Eq.  (35)  is 
used  in  Eq.  (37),  one  obtains 

fi\/  P-i  \  I  .2 


rE^11  ♦?’(»  -  ir)  h  ■  K 


-  k2)  1  +  (o2  -  k2)6  .  (39) 

m  mn  n  n  n,m  v  ' 


The  left-hand  side  of  Eq.  (39)  can  be  simplified  with  the  following  expression 
,S  1  *-*•  n  i  /  P; 
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mn 


iy h.  p.4.(i)  ^  (i  -  - l~\ 

2Z-»  i  V  P.  .  /  (z=H. 

i  11  1 


(40) 


The  function  is  the  symmetric  part  of  the  coupling  coefficient  defined 
in  Eq.  (6).  Equation  (40)  is  a  generalization  to  multiple  interfaces  of 
the  result  given  in  Appendix  A.  The  use  of  Eq.  (40)  in  Eq.  (39)  yields 
the  result,  • 

X  c  n  7  2  2 

2  —  RS  =  (a  -  k  1 1  ♦  (o  -  k  )S 

2  IT  mn  '  n  V  mn  v  n  nJ  n,m 
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For  m=n  Eq.  (41)  gives 

X 

O2  -  k2  =  -2  ^  B  +  0(H2)  , 

n  n  X  nm 

n 

since  I  <*0(11).  For  m^n  Eq.  (41)  gives 
inn 

-2X  BS 

j  -  n  mn 

mn  X  (k2  -  k2) 
n  n  m' 

2  2 

(In  obtaining  Eq.  (43),  °n  =  +  0(H)  was  employed.) 


(42) 


(43) 
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Now  return  to  Eq.  (32)  and  use  Eq.  (33)  to  get 

,  X 


E(*„  *  T’  *  <\)K  -  £[(*„  •  *  °l  K)  %  *  (ri, .  (♦, .  6*n) 

*  \(  \  *  6*„)]  • 

If  Eq.  (44)  is  multiplied  by  P<t>m  and  integrated  over  depth,  one  obtains 
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(45) 


In  obtaining  Eq.  (45),  terms  involving  6^  and  have  been  dropped  since 
•  2  >3 

they  are  of  order  H  and  H  ,  respectively.  Since  a  first  order  treatment 

of  the  field  is  desired,  Eq.  (45)  may  be  further  simplified  by  dropping 

•  2 

some  additional  terras  of  order  rl  and  higher.  To  do  this  note  that 
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n  n 


orn) 


i  =  0(in 

mn  1  - 


(46) 


With  Eq.  (46),  Eq.  (45)  becomes 

X 


*,  *  F1  *  k»X„  ■  -  °»>  X»  -E{(2Xn  *  r)  “»n  *  V„„}  ' 


If  Eq.  (42)  is  used  in  Eq.  (47),  the  result  is 

x  ,  ^  (\  .  N 

Y  +  _5L  +  k2  X  =  -  B  l  —  +  2X 
Xm  r  m  m  mn  \  r  ny 


where  X  /r  <<:  2X  has  been  assumed, 
n  n 


n^m 


(48) 


In  Eq.  (48)  the  coupling  coefficient  Amn  defined  by  Eq.  (7)  has  been 

dropped.  In  situations  in  which  the  only  range  dependence  of  the  medium 

•  2 

is  through  sloping  boundaries,  A  is  directly  proportional  to  H  and  should 

mn 

be  dropped  in  Eq.  (47).  If  the  medium  has  range  dependent  sound  speed, 

.  2  *2 

for  example,  A  will  have  terms  of  order  c  ,  cH,  and  H  .  In  this  case 
mn 
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•  •  •  2  *  * 
one  is  in  principle  entitled  to  retain  the  terms  of  order  c  and  cH  in 

Eq.  (48).  Since  values  of  c/c  that  one  might  expect  for  the  ocean  are 

of  the  same  size  or  less  than  typical  values  of  H/H  on  continental  slopes, 

it  is  asserted  that  Eq.  (48)  is  adequate  even  when  there  is  range  variation 

of  the  geoacoustic  parameters  of  the  waveguide. 

It  is  interesting  to  note  that  the  effect  of  requiring  the  boundary 
conditions  to  be  satisfied  to  order  tf  is  to  remove  the  diagonal  elements 
of  Bmn  (compare  with  Eq.  (2)).  Equation  (40)  reveals  that  is  completely 
off-diagonal  only  when  H=0  or  when  the  sloping  interfaces  are  pressure 
release  with  zssH^(r)]  =0 .  These  two  situations  correspond  to  cases 

in  which  conventional  coupled  mode  theory  is  capable  of  exactly  satisfying 
the  imposed  boundary  conditions.  The  occurrence  of  the  diagonal  elements 


of  B  is  therefore  directly  linked  to  the  use  of  improper  boundary  conditions 
at  sloping  interfaces. 


Return  now  to  the  equation  for  Equation  (42)  when  used  in  Eq.  (34) 


yields 


/a2  .2  2  \  r  , 

(  a  +  k  -  o  )  6<f) 

W  "/  n 


B  <p 
nn  n 


(49) 


Only  the  particular  solution  to  Eq.  (49)  is  desired  since  it  is  the  one 
that  vanishes  as  H  goes  to  zero. 

To  summarize,  a  coupled  mode  theory  expression  for  the  field  which 
satisfies  the  physically  correct  boundary  conditions  to  first  order  in 
the  bottom  slope  is  given  by  the  following 

Kr.z)  =£  Xn(r)  <yz;r)  ,  (50) 

n 

where  satisfies  Eq.  (48)  and  Eq.  (31).  Computation  of  the  field 
using  Eq.  (50)  as  opposed  to  Eq.  (1)  is  only  slightly  more  complicated. 
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Solving  for  X  is  no  more  involved  than  solving  for  the  R  of  conventional 
m  m 

coupled  mode  theory.  An  added  complication  arises  in  the  calculation  of 


♦  since  it  entails  the  computation  of  two  sets  of  functions  0  and  60  . 
n  n  n 

However,  the  additional  effort  to  calculate  the  60  may  not  be  as  great 

n 

as  it  first  appears.  An  approximation  to  Eq.  (49)  may  be  obtained  by  ignoring 

2  •  • 
the  first  order  part  of  o  and  using  a  WKB  approximation  for  X  (X =  ik  X  ) 

n  °  rr  nnnn 


to  obtain 


2 

9  2  *> 

— -  +  k  -  k‘  60  =  -  2  i k  B  0 

3Z*-  n  n  n  nn  n 


Once  the  0^  and  are  known,  the  <50^  may  be  obtained  by  a  routine  integration 

of  Eq.  (51).  The  only  eigenvalue  problems  to  be  solved  will  be  for  {k  , 

n 

0^}  at  whatever  range  values  are  required  in  the  solution  of  Eq.  (48). 

A  final  point  concerning  calculation  of  the  field  using  Eq.  (50)  is 
one  of  energy  conservation.  In  Appendix  B  it  is  shown  that  the  time  rate 
of  energy  flow  in  the  radial  direction  is  a  constant  to  first  order  in 
H  for  the  field  given  by  Eq.  (50).  Thus,  the  introduction  of  the  corrections 
proposed  here  not  only  result  in  a  theory  which  is  valid  to  a  specified 
order  in  H,  and  in  which  the  physically  correct  boundary  conditions  are 
satisfied,  but  also  in  a  theory  which  properly  conserves  energy  in  the 
case  of  zero  volume  absorption. 

IV .  NUMERICAL  CALCULATIONS 

In  this  section  some  numerical  calculations  illustrating  the  major 
points  of  this  paper  are  presented.  The  range  variable  waveguide  model 
chosen  for  these  calculations  is  the  isovelocity  wedge  shown  in  Fig.  2. 

The  surface  at  z=0  is  taken  as  pressure  release  and  the  sloping  bottom 
as  a  rigid  reflector.  For  this  waveguide  structure  the  normalized  depth 
functions  of  Eq.  (3)  are  given  by 


(52) 
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For  an  inicial  depth  of  Hq  *  500  m,  a  sound  speed  of  1.5  km/sec,  and  a 
frequency  of  20  Hzt  there  are  13  propagating  modes  present  at  the  source 
location. 


Since  the  mode  functions  and  coupling  coefficients  are  known  analytically, 
the  main  part  of  the  numerical  calculation  involves  computation  of  the 
radial  functions  of  Eq.  (48).  To. solve  Eq.  (48)  an  iterative  procedure 


was  employed  in  which  Eq.  (48)  is  replaced  with 

j(°) 

x(o)  *  Ji  ♦  k2  x<°)  -  n 
m  r  m  m  ~  0 


r*(i+l) 


♦  k2  X(i  +  1)  =  -Eb  -2—  +  2X(i) 


Equation  (53)  is  the  adiabatic  approximation*  to  the  set  of  coupled 


radial  equations  and  Eq.  (54)  is  an  approximation  to  Eq.  (48)  in  which 


the  right-hand  side  is  formed  from  the  previous  iterate.  The  iterative 

process  of  Eq.  (54)  is  continued  until  the  difference  between  two  successive 

iterations  is  less  than  some  specified  tolerance.  In  the  calculations 

presented  here,  convergence  to  four  significant  figures  was  required. 

The  iterative  procedure  just  described  is  very  similar  to  one  used  in  atomic 

14 

physics  applications  to  calculate  quantum  mechanical  bound  state  and 
scattering  state^  atomic  wave  functions. 

In  solving  for  the  radial  functions,  the  initial  conditions  for  the 
were  obtained  by  matching  to  the  known  horizontally  stratified  radial 
functions  at  r=rg  in  Fig.  2.  The  field  at  the  source  was  approximated 
as  having  only  outgoing  waves,  and  the  radial  functions  are  found  to  be 
proportional  to  Hankel  functions  of  the  first  kind  (for  e  lWt  time  dependence). 


X 

m 
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Figure  3  shows  the  result  of  a  computation  of  the  time  rate  of  change  of 

the  total  energy  in  the  radial  direction  as  a  function  of  range  using  Eq.  (B8). 

Also  shown  in  Fig.  3  is  the  same  calculation  performed  within  the  formalism 

of  conventional  coupled  mode  theory,  employing  the  vanishing  of  3i|V3z  at 

z=H(r).  The  stair-step  nature  of  the  top  curve  in  Fig.  3  is  a  consequence 

of  the  neglect  of  backscattered  energy  in  solving  for  the  radial  functions. 

The  abrupt  drops  in  the  top  curve  correspond  to  ranges  at  which  a  given 

mode  reaches  cutoff.  With  the  neglect  of  the  backscattered  field  on<.  is 

examining  only  the  time  rate  of  change  of  energy  flowing  away  from  the 

source  and  this  quantity  abruptly  changes  whenever  a  mode  drops  out.  The 

flat  portion  of  the  top  curve  in  Fig.  3  between  successive  mode  cutoff 

ranges  indicates  a  constant  rate  of  energy  transport  over  regions  having 

the  same  number  of  propagating  modes.  The  small  ripples  in  the  top  curve 

between  successive  mode  cutoff  ranges  arise  because  of  second  and  higher 

order  terms  that  are  retained  through  the  use  of  Eq.  (B8).  Although  the 

radial  functions  contain  only  zero  and  first  order  terms,  the  products 

•  2  -3 

of  radial  functions  in  Eq .  (B8)  give  rise  to  terms  of  order  H  and  H 
which  cannot  easily  be  removed.  The  dotted  curve  in  Fig.  3  is  a  conventional 
coupled  mode  theory  version  of  the  solid  curve.  The  nonconservation  of 
energy  effect  is  readily  apparent. 

Figure  4  shows  the  time  rate  of  change  of  energy  in  the  radial  direction 
associated  with  mode  No.  1.  This  was  obtained  from  Eq.  ( B8 )  by  eliminating 
the  summation  over  n  and  taking  n  to  be  1 .  In  Fig.  4  the  lower  curve  is 
a  conventional  coupled  mode  theory  evaluation  of  the  same  quantity.  Again, 
the  nonconservation  of  energy  is  easily  observed. 
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Figure  5  shows  a  comparison  between  the  radial  functions  obtained 
from  Eq.  (48)  and  the  adiabatic  approximation  to  Eq.  (48)  for  mode  No.  0. 
In  Fig.  5  the  1/^T  range  dependence  has  been  removed.  Note  the  changes 
in  amplitude  of  the  solid  curve  arising  from  energy  interchange  between 
modes,  a  process  not  present  in  adiabatic  mode  theory. 

Figure  6  shows  a  comparison  of  <j>n  and  6^^  for  n=*3.  The  curve  for 
was  obtained  from  the  approximation  of  Eq.  (51)  whose  solution  can 
be  found  analytically, 


The  range  at  which  the  computations  shown  in  Fig.  6  were  performed  was 

r=r  (see  Fig.  2)  which  was  taken  to  be  1.0  km. 
s 

Figure  7  depicts  the  absolute  square  of  the  velocity  potential,  which 
is  proportional  to  the  intensity,  computed  from  Eq.  (50)  for  the  waveguide 
structure  of  Fig.  2.  For  this  calculation  the  bottom  slope  was  5.0°  and 
the  source  and  receiver  depths  were  200  m  and  100  m,  respectively.  Note 
the  differences  between  the  adiabatic  and  coupled  mode  descriptions  of 
4*  ip.  Figure  8  shows  the  differences  between  the  two  curves  in  Fig.  7  plotted 
versus  range.  Though  large  differences  between  the  two  curves  exist,  the 
mean  difference  over  the  range  of  Fig.  8  is  +0.083.  This  leads  one  to 
suspect  that  if  a  range  averaged  description  of  acoustic  propagation  is 
desired,  then  mode-mode  coupling  effects  might  become  much  less  important 
or  even  insignificant. 

Figure  9  shows  the  phases  associated  with  the  two  velocity  potentials 
used  to  generate  Fig.  7.  It  should  be  noted  that  the  adiabatic  phase  accurately 


follows  the  phase  of  the  field  computed  including  mode-mode  coupling. 

The  approximate  agreement  of  the  two  phases  tends  to  support  the  conclusion 
of  a  previous  publication^  that  the  phase  structure  of  the  acoustic  field 
is  given  quite  accurately  by  the  adiabatic  approximation. 

V.  CONCLUSIONS 

Two  major  conclusions  may  be  drawn  from  this  work.  First,  the  conventional 

theory  of  sound  propagation  via  coupled  normal  modes  involves  approximations 

that  are  inconsistent  when  the  theory  is  applied  to  multiple  fluid  media 

having  range  variable  boundaries.  A  consequence  of  this  inconsistency 

is  nonconservation  of  energy.  Second,  a  correction  to  the  conventional 

theory  can  be  derived  such  that  the  field  so  obtained  satisfies  the  physically 

proper  boundary  conditions  to  first  order  in  the  slope  of  range  variable 

boundaries  that  are  present.  With  these  first  order  corrections,  energy 

also  becomes  conserved  to  first  order.  For  a  more  detailed  account  of 

some  aspects  of  the  work  presented  in  this  article,  see  Ref.  17. 

The  most  important  difference  between  the  conventional  theory  and 

the  corrected  theory  is  the  additional  quantity  6$^  that  must  be  computed. 

Depending  on  the  type  of  application,  the  added  computational  difficulty 

imposed  by  the  6^  may  be  negligible  compared  to  the  other  calculations 

that  must  be  performed,  e.g.,  solving  Eq.  (48).  The  are  not  needed 

until  the  field  is  to  be  constructed  from  the  three  components  X  ,  <p  , 

n  n 

and  5(f)  .  This  is  not  true  of  the  4>  which  must  be  known  as  a  function 
n  n 

of  range  so  that  the  coupling  coefficients  may  be  determined.  Therefore, 
if  one  is  interested  in  the  field  as  a  function  of  depth  at  one  particular 
range,  then  only  the  <5$  corresponding  to  that  range  need  be  calculated 
with  a  minimal  increase  in  computational  effort. 
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If  the  field  as  a  function  of  range  is  desired,  then  the  computation 
of  the  can  become  more  costly  in  terras  of  effort  depending  on  the  range 
grid  on  which  the  field  is  required  to  be  computed.  One  particular  type 
of  calculation  using  the  corrected  theory  that  requires  no  extra  computational 
effort  arises  in  the  computation  of  the  random  phase  propagation  loss  defined 

by  i 1 2 


i*r 


rp  ~  y  x  x  $  $ 

n  n  n  n 


Since  by  approximation  Eq .  (51)  6(J>n  is  pure  imaginary,  Eq.  (56)  to  fi-st 


order  becomes 


kl2on~  y  x*  x  <p  <$> 

RP  *-•  n  n  n  Tn 


and  the  6<J>n  need  not  be  computed  at  all.  This  could  be  very  significant 
for  a  variety  of  practical  applications  where  only  an  estimate  of  average 
field  level  is  desired. 
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APPENDIX  A 


This  appendix  is  devoted  to  the  derivation  of  several  equations  appearing 
in  Section  II. 

The  rate  of  change  of  power  radiated  through  the  surface  of  a  semi¬ 
infinite  cylinder  was  found  to  be,  (Eq.  (15)), 


Upon  introducing  the  partially  separated  solution,  Eq .  (1),  one  has 


or, 

3P 

37 


If  =  -mini  [pj'PjH  (i>\  -  ^)]  •  (A1) 

Upon  introducing  the  partially  separated  solution,  Eq .  (1),  one  has 

(A2) 

or, 

=  j  iiruir  53  [p  H  (r  R*  -  R  R*\ 
dr  2  *-•  1  V  n  m  m  n  /  V  m  m  / 

m,n  L  '  ' 

+  p  h(r  R*  -  R*  R  ^  ^ ^ 1  ^  (V1)  -  )  (A3) 

1  \  n  m  n  m )  n  \  m  m  /J  . 

Now,  in  the  second  term  of  Eq.  (A3)  we  introduce  the  notation  for  the 

quantity  R  R*  -  R*R  ,  and  in  the  first  term,  after  interchanging  the  summation 

indices  (n,  ml,  one  finds, 

&  »  i  imr  £  fp,  H  (r  R*  -  R  R* )  -  tm) 

9r  4  \nm  m  n  J  Yn  Vn  ym  / 

+  p  H  (r  R*  -  R  R*\  (<J>(1)  -  <f>(2))l 

1  \mn  nm/  m  \rn  n  /J 

+  \  imur  53  [p,  A  Q  -  4>(2'l)l  •  (A4) 

2  jffo  11  n  \rm  m  / 

Upon  using  the  continuity  relation  Pj<t>^^=  the  first  term  in  Eq.  (A4) 


can  be  rewritten  to  obtain 


5?  ■  r  i"“r  £„{[pi  ”  *il)  (‘  -  DiA)]  p„.} 

*  5  {  [°>  "  ^  W1’  -  .  (AS) 
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where  we  have  introduced  the  notation  P  for  the  combination  of  radial 

non 

functions  appearing  in  Eq.  (A5)  (see  Section  II).  It  now  remains  to  express 
the  two  factors  in  brackets,  [•],  in  terms  of  mode  coupling  coefficients. 

The  depth  functions  4>n(z)  have  been  assumed  to  form  a  complete  orthonormal 
set  obeying  ^ ( 0 )=0 ,  continuity  of  P4>r  and  ^  at  the  interface  z*»H(0, 
and  a  radiation  condition  as  Z-*00.  The  orthogonality  relation  for  these 
functions  can  be  differentiated  with  respect  to  r  to  obtain 

fr  H(r)  oo  i  i 

Pi  /  dz  +  P2  /  dz  f  =  IF  6nm  =  °  '  (A6) 

o  H(r)  J 


P1  H  ^  ^m13  -  p2  H  ^  +  /  dz  p(^n  K  +  K  =  0  *  (A7) 

o 

Upon  introducing  the  definition  of  the  B  into  Eq.  (A7),  and  using  the 
continuity  of  p<f>n>  we  have 


p.  H  <)>(1)  fl  -  P,/6~)  +  B  +  B  =  0. 

1  n  rm  \  V  2)  nm  mn 


The  desired  result,  Eq.  (20),  is  thus  obtained  at  once. 


'  T  C1  (J  -  »A)  • 


It  should  be  noted  that  the  symmetric  part  of  B  is  proportional  to  H. 

nm 

Thus,  in  cases  such  as  those  originally  considered  in  Refs.  1  and  2,  B 
is  indeed  antisymmetric  as  claimed.  The  symmetric  part  of  Bmn  is  also  proportional 
to  the  density  discontinuity;  hence,  it  might  therefore  be  anticipated 
that  for  small  slopes  and  high  porosity  sediments,  such  as  clays  and  silts, 

S 

B  would  be  small, 
nm 


The  corresponding  expression  foi  the  second  term  in  Eq.  (A5)  can  be 

obtained  by  differentiating  the  defining  equation  for  B  (r), 

nm 

3Bnm  9  f  9<^m  f  9<f>m 

9r  ~  3r  P1  J  dz  TT  +  p2 /dz  K  IF  ' 

o  H 


which  yields 


Bnm  =  P1  H  ^  *m'>  '  p2  H  ^  +  /  dz  p(*n  +  *n  $m)  »  (/ 


or,  after  introducing  the  definition  of  A  t 

nm 


*  *1"  (V  ■  C)  */-  p  K  ♦.  ^ 


Thus,  the  matrix  B  -A  is  given  by, 
nm  nm 


B  -  A 
nm  nm 


-  p.  H  ♦  <*)  (*»>  .  *C2J\  (  .  . 

1  "  V^m  %  )  J  dz  p  w 


It  is  easily  seen  that  the  second  term  on  the  right-hand  side  of  Eq.  (A13) 
is  completely  symmetric  in  (m,n),  and  consequently  the  desired  antisymmetric 
matrix  is  given  by 


(5n.  -  vf  *0,  «  [*?>  (if>  -  »«)  -  ♦  «>  (}»>  -  }<2>)  . 

^  '  - 


which  is  the  result  previously  quoted  in  Eq.  (19). 


(A14) 


APPENDIX  B 


In  this  appendix  it  is  shown  that  the  rate  of  radial  energy  flow  for 
the  field  given  by  Eq.  (50)  is  a  constant  to  first  order  in  the  bottom 
slope  H. 


The  time  averaged  radial  energy  flux  is  given  by 

Jr  =  I  Re[p  vr]  ' 


(Bl) 


where  p  is  the  acoustic  pressure  and  v^  the  radial  particle  velocity  defined 


by 


p  =  iupip 


dtp 

V  =  TT  . 

t  dr 


(B2) 


In  the  following  developments  it  is  useful  to  remove  the  1  /v/r~r_  dependence 
from  the  radial  functions.  In  this  appendix  the  field  will  be  assumed 
to  have  the  form 

G  (r) 

(z:r)  • 


(B3) 


where  the  equation  for  is  obtained  from  Eq.  (48). 

‘  (k"  *1?)  c"  ’  -?  2B»»  6 
m/n 


(B4) 


with  G  (r)=\Zr~  X  (r). 
n  n 

An  expression  for  the  time  averaged  radial  energy  flux  may  be  obtained 
from  Eq.  (Bl)  through  Eq.  (B3). 


* 

J,  -  f  Im<!  E  k  s*  -  ♦ „.c*  6  ♦  i* 

1  m,n  \  nm  iT  I  m  n  mnnm 


(B5) 
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The  time  rate  of  energy  flow  in  the  radial  direction  is  obtained  from  Eq. 

by  integrating  over  a  cylindrical  surface  of  radius  r.  This  procedure 

yields  the  radial  rate  of  energy  flow,  P. 

00 

P  =  J  2vt  dz 
0 


P  =  -m  Im  J  £  (g  G*  -  Gm  G 


(«  +  I*  +  I  ) 

2r  J  \  n,m  nm  mn/ 


+  G*  G  B 
m  n  nm 


(B5) 


(B6) 


In  Eq.  (B6),  I  is  given  by  Eq.  (43)  and  terms  of  higher  order  than  H 
mn 

were  neglected.  If  one  uses  the  approximation  Xfc*(ik  )X 

n  n  n 

in  Eq.  (43),  1*  +1  is  seen  to  be  pure  imaginary.  Therefore,  Eq.  (B6) 
nm  mn 

becomes 


P  = 


G 

n 


G 


m 


+ 


★ 


G 

n 


+  £  (g  G*  -  G*  G  )  B  1 
m7n  '  n  m  n  m  /  nm  f 


(B7) 


In  Eq.  (B7)  the  terms  G*G  /2r  were  neglected  in  comparison  to  G  G*.  If 


m  n 


n  m 


G  «ik  G  +0(fi)  is  used  in  the  second  and  third  terms  of  Eq.  (B7),  the  result 
n  n  n 


is 


P  = 


{e(g  G  -  G*  6  )  *£(g  G*  -  G*  G  )  rB  -  ik  (i*  +  I  n)ll 
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(B8) 


In  Eqs.  (B6)  and  (B7)  the  approximation  X^= ik^X^  was  employed.  This 

approximation  gives  an  expression  for  Xn  that  is  correct  through  terms 

of  zero  order  in  ft  and  is  appropriate  whenever  terms  of  first  order  in 

H,  such  as  B  ,  occur  as  multiplicative  factors  of  ft  •  This  approximation 
mn  r  n 

may  be  viewed  as  a  WKB  expression  for  the  derivative  of  the  forward  going 
component  of  the  adiabatic  solution  of  the  radial  equation.  The  backscattered 
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component  can  be  shown  to  be  of  order  H  compared  to  the  forward  going  component 
and  is  ignored. 


To  show  that  P  is  independent  of  range,  Eq.  (B8)  is  differentiated 
with  respect  to  range.  This  yields 


p‘  '  (‘„  V  <  *„) 
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ik  (i*  +  r 

m  \  nm  mn> 


)]} 


(B9) 


In  obtaining  Eq.  (B9)  terms  of  order  H  were  not  differentiated  since  they 
would  produce  quantities  of  order  H  .  From  Eq.  (B4)  it  follows  that 

~  (g  G*  -  G*  G  )  =  £  2B  (g*  G  -  G  G*)  (BIO) 

dr\nn  nn/„nm\nm  nm/ 

Equation  (BIO),  (B9)  gives 
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If  one  uses  the  WKB  approximation  in  Eq.  (43),  I*  +  I  is  seen  to  be 

rr  ^  nm  mn 


(Bll) 
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(B12) 


The  use  of  Eq.  (B12)  in  Eq.  (Bll)  yields 
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=  {  £  i(k  -  k  )  BS  (g  G*  +  G*  G  \  1 

r  2  [  njm  '  n  ">/  mn  \  n  m  n  m/  f 
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(B13) 


Equation  (B13)  vanishes  by  virtue  of  the  symmetry  properties  of  the  quantity 
being  summed.  Therefore,  P  is  zero  through  terms  of  order  H, 

Jr  =  0  ♦  0 (h 2)  .  (B1 4 ) 
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FIGURE  2 

ISOVELOCITY  WEDGE  SHAPED  WAVEGUIDE 
ASSUMED  FOR  NUMERICAL  COMPUTATIONS 
c  =  1500  m/sec,  Hq  =  500  m 


ARL:UT 
AS-79-1950-P 
SRR  -  GA 
9-  13  -79 
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TOTAL  POWER 


FREQUENCY:  20  Hz 
BOTTOM  SLOPE:  2.5° 

WITH  FIRST  ORDER  CORRECTIONS 


FIGURE  3 

TOTAL  POWER  TRANSPORTED  IN  RADIAL 
DIRECTION  (W)  versus  RANGE 

THE  ARROWS  DENOTE  MODE  CUTOFF  RANGES 
SOURCE  STRENGTH  =  1  m3/sec,  SOURCE  DEPTH  =  200  m 


ARL:UT 
AS-79-1951-P 
SRR  * GA 
9-13-79 
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FIGURE  4 

POWER  TRANSPORTED  BY  MODE  1  (W)  versos  RANGE 
SOURCE  STRENGTH  =  1  m3/sec,  SOURCE  DEPTH  =  200  m 


FREQUENCY:  20  Hz 
BOTTOM  SLOPE:  2.5° 

-  WITH  FIRST  ORDER  CORRECTIONS 

.  ADIABATIC  APPROXIMATION 


n 

10 
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RANGE  -  km 


FIGURE  5 

REDUCED  RADIAL  FUNCTION  (V^X~n) 
FOR  MODE  0  versus  RANGE 

SOURCE  STRENGTH  1  m3/sec, 
SOURCE  DEPTH  200  m 


-  ADIABATIC  MODE  THEORY 

- COUPLED  MODE  THEORY  WITH  CORRECTIONS 


RANGE  -  km 


FIGURE  7 

ABSOLUTE  SQUARE  OF  VELOCITY  POTENTIAL  versus  RANGE 
SOURCE  STRENGTH  -  1  m3/sec,  SOURCE  DEPTH  -  200  m 


DIFFERENCE  IN  4>*4>  -  km4/sec2 
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i  The  effect  of  sediment  rigidity  on  bottom  reflection  loss  in  a  typical 

deep  sea  sediment 

Paul  J.  Vidraar,  Applied  Research  Laboratories,  The  University  of  Texas 

at  Austin,  Austin,  TX  78712 

(Received 


The  effect  of  sediment  rigidity  on  bottom  reflection  loss,  RL  *  -20  log^l  r| 
where  R  is  the  plane  wave  reflection  coefficient,  is  studied  using  a  com¬ 
putational  model.  A  single  inhomogeneous  turbidite  layer  overlying  a  homogeneous 
substrate  is  treated  with  an  emphasis  on  low  frequencies  (10  to  200  Hz) 
and  low  grazing  angles  (0°to  45°).  We  find  that  sediment  rigidity  can 
be  neglected  for  thick  layers  (-500  m)  while  it  can  produce  significant 
increases  in  RL  for  thin  layers  (-36  m) .  The  frequency  dependence  of  RL 
for  thin  layers  has  a  low  frequency  regime  characterized  by  large  peaks 
(-25  dB)  recurring  at  short  (-2  Hz)  intervals.  At  high  frequencies  there 
remains  a  persistent  increase  in  RL  (-4  dB).  The  most  important  mechanism 
for  sediment  S  wave  excitation  is  found  to  be  compressional  wave  conversion 
at  the  sediment-substrate  interface. 

PACS  numbers:  43.20.Fn  43. 30. Dr 
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INTRODUCTION 


The  acoustic  interaction  with  the  ocean  bottom  is  an  important  factor 

in  determining  the  propagation  of  low  frequency  sound  waves  in  the  oceans. 

Because  of  the  long  wavelengths  involved,  the  acoustic  energy  penetrates 

into  the  ocean  bottom  and  interacts  with  the  subbottom  structure.  The 

compressional  (P),  shear  (S),  and  interface  waves  excited  in  the  sediment- 

substrate  structure  are  important  energy  loss  mechanisms  influencing  low 

frequency  sound  propagation  in  the  oceans.  A  useful  and  traditional  measure 

of  this  bottom  interaction  is  the  plane  wave  reflection  coefficient  R  and 

its  analog  the  reflection  loss,  RL  *  -20  log^l  r|  . 

Theoretical  studies  for  an  assumed  fluid  sediment,  using  numerical 

models  to  compute  R,  have  investigated  the  influence  of  several  subbottora 

features  on  RL  at  low  frequencies.  Density  gradients*  and  the  roughness 

2 

of  the  water -sediment  interface  have  a  negligible  effect  for  typical  deep 

3  4  5 

sea  sediment  types.  The  effect  of  P  wave  speed  ’  and  absorption  profiles 
have  also  been  treated.  Stoneley  wave  excitation  at  the  interface  between 
the  (fluid)  sediment  and  a  (solid)  substrate  has  been  identified  as  a  loss 
mechanism.*’’  ^ 

The  effect  of  sediment  rigidity  (shear  wave  propagation)  on  R  has 

received  little  attention  and  is  not  well  understood.  Although  included 

8  9 

in  some  numerical  models,  ’  no  systematic  studies  have  been  made  to  deter¬ 
mine  the  importance  of  energy  lost  to  sediment  S  waves.  A  detailed  study 
at  this  time  would  be  useful  since  parameter  ranges  can  now  be  constrained 
by  estimates  of  sediment  S  wave  properties  that  have  recently  become  avail¬ 
able.*^  *^  Such  a  study  is  important  because  the  relatively  high  attenuation 

of  sediment  S  waves  make  their  excitation  a  potentially  important  energy 

.  it  14  .  .  . 

loss  mechanism.  Recent  work  ’  indicates  that  sediment  S  wave  excitation 

can  be  the  dominant  energy  loss  mechanism  is  thin  sediment  layers. 
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In  this  paper  we  investigate  the  effect  of  sediment  rigidity  on  R 
for  a  typical  deep  sea  subbottom  structure  composed  of  a  single  inhomogeneous 
turbidite  layer  overlying  a  homogeneous  basalt  substrate.  Our  emphasis 
is  on,  but  not  restricted  to,  low  grazing  angles  and  low  frequencies. 

The  effect  of  sediment  rigidity  is  determined  by  comparing  RL  for  a  solid 
sediment  and  a  fluid  sediment.  The  difference  between  them,  ARL,  i8  due 
to  sediment  rigidity.  We  find  that  the  dependence  of  Arl  on  the  thickness 
of  the  sediment  layer  leads  to  an  empirical  classification  of  sediments 
into  "thick"  and  "thin"  categories.  For  a  thick  (518  m)  layer  Arl=q  at 
low  grazing  angles  and  sediment  rigidity  can  be  neglected.  For  a  thin 
(36  ra)  layer  at  a  grazing  angle  of  20°,  ARL  ~ 20  dB  at  20  Hz  and  ARL?4  dB 
at  200  Hz;  sediment  rigidity  in  this  case  is  not  negligible  but  produces 
significant  increases  in  RL.  The  frequency  dependence  of  RL  for  the  36  m 
thick  layer  reveals  high  and  low  frequency  regimes  with  strikingly  different 
behavior.  For  a  grazing  angle  of  20°,  ARL  is  about  2  dB  at  low  frequencies 
with  very  large  peaks  (~25  dB)  occurring  at  2  Hz  intervals.  The  peaks 
decrease  in  magnitude  at  higher  frequencies  until  they  entirely  disappear. 

At  high  frequencies  there  remains  a  persistent  4  dB  increase  in  RL  compared 
to  that  computed  for  a  fluid  sediment.  Our  study  suggests  that  the  mechanism 
for  sediment  S  wave  excitation  is  primarily  P  wave  conversion  at  the  sediment- 
substrate  interface.  The  gradient-driven  P-S  conversion  is  negligible 
for  these  sediment  types  for  frequencies  above  about  3  Hz.^  This  excitation 
mechanism  explains  the  classification  of  sediment  layers  into  thick  and 
thin  categories.  For  a  thick  sediment  the  P  wave  has  its  turning  point 
well  above  the  sediment-substrate  interface.  Consequently,  sediment  S 
wave  excitation  is  negligible  and  the  sediment  can  be  accurately  treated 
as  a  fluid.  For  a  thin  sediment  the  P  wave  significantly  interacts  with 
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the  substrate.  Shear  waves  are  excited  in  the  sediment,  and  the  sediment 
must  be  treated  as  a  solid. 

The  separation  into  high  and  low  frequency  regimes  in  thin  layers 
is  due  to  interference  effects  related  to  S  wave  propagation  within  the 
sediment.  At  low  frequencies  sediment  shear  wave  attenuation  is  small 
and  interference  effects  will  dominate  the  dependence  of  R  on  frequency. 

Large  peaks  in  RL  will  occur  with  a  separation  in  frequency  related  to 
changes  in  S  wave  phase  of  2tt.  At  higher  frequencies  the  increased  shear 
wave  attenuation  decreases  the  magnitude  of  this  peak  structure.  At  high 
frequencies  S  waves  excited  at  the  sediment-substrate  interface  are  totally 
absorbed  within  the  sediment,  resulting  in  an  almost  constant  additional 
loss  compared  to  the  RL  of  a  fluid  sediment. 

The  tool  used  in  this  study  is  a  recently  developed  computational 
model^  of  R  which  allows  sediment  shear  wave  propagation.  The  model  assumes 
horizontal  stratification  and  treats  a  single  sediment  layer  overlying 
a  seraiinf inite  homogeneous  substrate.  The  model  is  based  on  numerical 
integration  of  the  depth-separated  wave  equations  for  the  potentials  giving 
rise  to  the  compressional  and  shear  waves  in  the  sediment  and  thus  includes 
all  wave  properties  such  as  penetration  beyond  turning  points,  boundary 
wave  excitation,  etc.  The  use  of  numerical  integration  permits  arbitrary 
depth  variations  of  sediment  parameters  and  allows  them  to  be  individually 
varied  to  determine  their  influence  on  R.  The  basic  approximation  in  the 
model  is  the  use  of  the  Helmholtz  equations  with  depth  dependent  wave  number 
to  describe  the  potentials.  This  is  esentially  a  high  frequency  approxi¬ 
mation  in  which  the  continuous  coupling  between  shear  and  compressional 
waves  and  effects  directly  dependent  on  gradients  are  neglected.  For  parameters 
typical  of  deep  sea  sediments,  this  model  is  accurate  for  frequencies  above 
about  10  Hz. 
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II.  REFLECTION  LOSS  OF  A  HYPOTHETICAL  TURBIDITE  LAYER 

In  this  section  we  consider  the  effect  of  sediment  rigidity  on  the 
reflection  loss  from  a  hypothetical  turbidite  layer.  RL  is  studied  as 
a  function  of  grazing  angle  ®,  frequency  f,  and  sediment  layer  thickness 
H  (the  grazing  angle  is  the  complement  of  the  angle  of  incidence  Y ,  i.e., 

0  =90-  y) . 

The  depth  structure  of  the  turbidite  layer  is  given  in  Table  I.  These 

.  9 

parameters  were  obtained  from  the  recent  work  of  Fryer.  Constant  gradients 
are  assumed  between  the  depths  given  in  Table  I.  The  parameter  values 
and  gradients  of  this  depth  structure  are  representative  deep  sea  turbidite 
layers.  The  attenuation  is  assumed  to  depend  linearly  on  frequency.^ 

For  our  computational  model  it  was  necessary  to  truncate  Fryer's  650  m 
sediment  layer  at  518  ra  to  produce  accurate  results  at  low  grazing  angles 
for  a  frequency  of  20  Hz.  For  a  significantly  thicker  layer  (650  m)  the 
large  S  wave  attenuation  combines  with  the  finite  computer  word  length 
to  produce  a  loss  of  precision  for  grazing  angles  below  about  20°. 

Figure  1  shows  RL  as  a  function  of  0  for  f  =  20  Hz  and  H  =  518  m. 

The  three  curves  were  obtained  for  different  sediment-substrate  configurations. 
The  dotted  line  was  obtained  for  a  fluid  sediment  and  fluid  substrate  (the 
FF  case),  the  dashed  line  for  a  fluid  sediment  and  solid  substrate  (the 
FS  case)  and  the  solid  line  for  a  solid  sediment  and  solid  substrate  (the 
SS  case).  The  remaining  case  of  a  solid  sediment  and  a  fluid  substrate 
(the  SF  case)  is  not  shown  in  Fig.  1.  It  is  essentially  identical  to  the 
FF  case  over  the  entire  range  of  grazing  angles,  frequencies,  and  sediment 
thicknesses  considered.  The  parameters  for  the  fluid  sediment  and  substrate 
were  obtained  from  Table  I  by  ignoring  the  S  wave  parameters. 
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A  comparison  of  Che  curves  of  Fig.  1  shows  Che  effecC  of  sediraenC 

rigidiCy  in  Chis  Chick  layer.  Case  FF  has  Che  expecCed  small  RL  for  small 

9  followed  by  an  oscillaCory  sCrucCure  due  Co  inCerference  beCween  Che 

wave  direcCly  reflecCed  from  Che  waCer-sedimenC  inCerface  and  Che  P  wave 

refracCed  in  Che  sedimenC  and  reCurning  Co  Che  waCer.  Beyond  Che  criCical 

angle  for  subsCraCe  P  wave  propagaCion  (8^=70*)  RL  increases  as  energy 

is  carried  away  by  subsCraCe  P  waves.  A  comparison  of  Che  FS  and  FF  cases 

shows  Che  effecC  of  subsCraCe  rigidicy  alone.  AC  boCh  small  and  large 

grazing  angles  Che  FS  and  FF  cases  are  nearly  idenCical  showing  ChaC  subsCraCe 

rigidicy  has  a  negligible  effecC  ac  Chese  angles.  RL,  however,  increases 

beCween  Che  criCical  angle  for  subsCraCe  S  wave  propagaCion,  6g,  and  6^. 

SubsCraCe  S  waves  carry  away  addicional  energy  in  Chis  angular  range. 

Comparing  Che  SS  and  FS  cases  shows  Che  addicional  effecCs  of  sedimenC 

rigidicy.  Excepc  for  grazing  angles  beCween  6  and  6  ,  sedimenC  S  wave 

s  p 

propagaCion  has  a  negligible  effecC  on  RL.  The  major  effecC  beCween  Che 

criCical  angles  is  a  shifc  in  Che  peak  sCrucCure. 

q  . 

Comparison  of  Fig.  1  wich  Fryer's  Fig.  8  shows  qualicacive  agreemenC 
between  the  corresponding  FF  and  SS  cases.  Fryer  obtained  his  result  using 
a  different  computational  technique  and  for  a  thicker  (650  m)  sediment 
layer.  Our  computation  of  the  intermediate  FS  case  allows  us  to  draw  the 
conclusion  that  for  this  particular  subbottom  structure  sediment  rigidity 
is  negligible  at  low  grazing  angles  and  for  a  frequency  of  20  Hz. 

Figure  2  shows  RL  for  a  thinner  sediment  layer,  H  =  252  m,  at  f  =  20  Hz. 
The  FS  and  SS  cases  are  shown.  Comparison  of  the  solid  (SS)  and  dashed 
(FS)  curves  shows  that  the  presence  of  sediment  rigidity  produces  the  expected 
shift  in  the  oscillatory  structure  between  the  critical  angles  9g  and  8^. 

In  addition,  the  SS  case  has  an  increased  RL  beginning  at  a  minimum  grazing 
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angle,  0q=240  ,  and  extending  up  to  about  45°.  This  increase  in  RL  is  entirely 
due  to  sediment  rigidity.  For  grazing  angles  below  0q  the  effect  of  sediment 
rigidity  is  negligible. 

Figure  3  shows  RL  for  the  252  m  layer  as  a  function  of  frequency  at 
0=2Oo<0o  for  the  frequency  range  from  10  to  100  Hz.  A  comparison  of  the 
SS  and  FS  cases  shows  that  for  frequencies  above  15  Hz  sediment  rigidity 
is  negligible.  Below  15  Hz  the  difference  between  the  curves  may  not  be 
meaningful  since  this  is  near  the  low  frequency  limit  of  the  validity  of 
our  computational  model. 

A  consideration  of  RL  as  a  function  of  0  for  sediment  thicknesses 

less  than  252  m  reveals  a  relationship  between  H  and  0  .  As  H  decreases, 

r  o 

also  decreases.  Thus,  as  sediment  thickness  decreases,  sediment  rigidity 

causes  changes  in  RL  at  steadily  lower  grazing  angles.  In  fact,  sediment 

S  wave  excitation  can  be  the  dominant  energy  loss  mechanism  in  thin  sediment 

layers.  This  is  illustrated  in  Fig.  4,  which  compares  RL  for  the  SS  and 

FS  case  at  H  =  36  m  and  f  =  20  Hz.  The  angle  0  is  0° .  Sediment  S  wave 

o 

excitation  is  responsible  for  the  dramatic  20  dB  increase  in  RL  near  0=12°. 

These  results  lead  to  an  empirical  classification  of  sediments  into 
thick  and  thin  layers  based  on  the  effect  of  sediment  rigidity.  For  thick 
sediments  little  energy  is  transferred  to  sediment  S  waves  and  the  sediment 
can  accurately  be  treated  as  a  fluid  at  low  grazing  angles.  For  a  thin 
sediment  significant  energy  can  appear  in  sediment  S  waves.  Sediment  rigidity 
is  not  negligible  and  the  sediment  must  be  treated  as  a  solid. 

The  dependence  of  RL  on  frequency  for  the  same  36  m  layer  is  shown 
in  Fig.  5.  The  SS  and  FS  cases  are  compared  from  10  to  200  Hz  at  0=20°  . 

The  SS  case  has  high  and  low  frequency  regimes  with  strikingly  different 
properties.  At  low  frequencies  the  effects  of  sediment  rigidity  increase 


90 


the  level  of  RL,  relative  to  the  FS  case,  by  about  3  dB  with  additional 
very  large  peaks  of  up  to  25  dB.  The  peaks  repeat  with  an  almost  constant 
separation  of  2.65  Hz.  The  magnitude  of  the  peaks  generally  decrease  with 
increasing  frequency;  however,  any  individual  peak  may  be  higher  or  lower 
than  neighboring  peaks.  The  peaks  disappear  by  about  90  Hz  and  the  high 
frequency  regime,  characterized  by  rather  broad  oscillations  in  RL,  begins. 
The  maxima  of  these  broad  oscillations  are  separated  in  frequency  by  about 
65.5  Hz  and  occur  in  both  the  SS  and  FS  cases.  In  this  high  frequency 
regime,  sediment  S  wave  excitation  produces  an  almost  constant  increase 
in  RL  of  about  4  dB. 

The  large  peaks  in  the  low  frequency  regime  of  Fig.  5  are  associated 
with  interference  effects  due  to  sediment  S  wave  propagation.  This  iden¬ 
tification  is  made  by  noting  that  the  frequency  interval  between  peaks 
is  related  to  S  wave  parameters.  One  expects  interference  effects  to  recur 
when  the  relative  phase  of  the  waves  involved  change  by  2 tt.  The  phase 
of  the  S  wave,  accumulated  in  traveling  through  the  sediment  and  back, 
changes  by  2 tt  for  a  frequency  change  Af=v,/2H,  where  v^  is  the  average 
sediment  S  wave  speed  (the  sediment  shear  speed  is  so  small  that  the  S 
wave  propagates  nearly  normal  to  the  plane  of  stratification).  The  measured 
separation  between  peaks  is  within  5%  of  the  estimated  value  of  Af.  The 
corresponding  quantity  for  P  waves  (including  the  grazing  angle  correction) 
estimates  a  frequency  change  of  65.4  Hz  which  is  a  factor  of  30  too  large 
to  explain  the  peak  separation.  In  addition,  the  decrease  in  peak  height 
with  frequency  appears  to  be  related  to  the  increased  attenuation  of  S 
waves  with  frequency.  At  90  Hz  the  attenuation  is  estimated  to  be  23  dB 
for  an  S  wave  traveling  one  way  through  the  sediment.  This  suggests  that 
S  wave  interference  would  be  negligible  above  90  Hz  as  observed  in  Fig.  5. 


The  gentle  oscillations  above  90  Hz  in  Fig.  5  are  due  to  P  wave  inter¬ 
ference  effects.  This  conclusion  is  reached  because  the  same  oscillations 


i 

are  observed  in  both  the  SS  and  FS  cases.  The  measured  separation  between 
maxima  is  65.5  Hz  which  agrees  very  well  with  the  estimated  Af  =  65.4  Hz. 

Figure  6  shows  that  sediment  rigidity  is  still  significant  in  the 
high  frequency  regime.  The  SS  and  FS  cases  are  compared  at  H  =  36  m  and 
f  =  200  Hz.  The  increase  in  RL  due  to  sediment  S  wave  excitation  is  signifi¬ 
cant  (~4  dB)  between  15°  and  30°.  The  minimum  angle,  6  ,  is  about  5°. 

Inspection  of  Fig.  5  near  20  Hz  shows  that  the  very  large  increase 
in  RL  seen  in  the  SS  case  in  Fig.  4  is  due  to  the  fortuitous  location  of 
one  of  the  interference  peaks  at  20  Hz.  For  comparison  with  Fig.  4,  RL 
at  21  Hz  (a  minimum  in  RL)  is  given  in  Fig.  7.  The  increase  in  RL  due 
to  sediment  S  wave  excitation  is  still  present  at  low  angles  but  its  magnitude 
is  2  dB  rather  than  20  dB.  This  is  still  an  important  increase  in  RL  when 
compared  to  the  RL=0  in  the  fluid  sediment  case. 

The  dependence  of  9q  on  H  gives  an  indication  of  the  physical  mechanism 
responsible  for  the  increase  in  RL  due  to  sediment  rigidity.  For  a  given 
grazing  angle  the  compress ional  wave  in  the  sediment  has  a  turning  point 
at  a  depth  Ht.  For  a  thick  sediment  layer  H  will  be  well  above  the  sediment- 
substrate  interface.  As  the  grazing  angle  increases,  H  increases.  At 
some  critical  angle  0  ,  Ht~H  and  the  compress ional  wave  interacts  with 
the  sediment-substrate  interface.  As  the  thickness  of  the  sediment  layer 
decreases,  the  angle  0^  also  decreases.  This  is  the  same  qualitative  behavior 
exhibited  by  the  angle  9  ,  above  which  sediment  rigidity  causes  an  increase 
in  RL.  This  similarity  suggests  that  sediment  shear  wave  excitation  by 
compress ional  wave  conversion  at  the  sediment-substrate  interface  produces 
the  increase  in  RL. 
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Two  other  observations  support  this  identification  of  P  wave  conversion 


at  the  sediment-substrate  interface  as  the  dominant  mechanism  for  sediment 
S  wave  excitation.  One  concerns  the  negligible  influence  of  the  water- 
sediment  interface  while  the  other  shows  the  importance  of  the  substrate 
interface.  S  wave  excitation  at  the  water-sediment  interface  appears  to 
be  negligible.  The  properties  of  the  water  interface  remained  constant 
while  the  sediment  thickness  varied  and  large  changes  in  RL  due  to  sediment 
rigidity  occurred.  This  suggests  the  lack  of  a  casual  relationship  between 
the  constant  water-sediment  interface  conditions  and  the  significantly 
increased  RL.  The  importance  of  the  substrate  interface  is  further  indicated 
by  the  negligible  difference  between  RL  for  the  FF  and  SF  cases  and  the 
large  changes  seen  between  the  SS  and  FS  cases.  This  indicates  that  the 
nature  of  the  substrate  is  important.  Since  the  substrate  parameters  enter 
our  model  only  through  the  sediment-substrate  interface  conditions,  this 
also  indicates  that  sediment  S  waves  originate  at  the  substrate  interface. 
III.  SUMMARY 

A  study  of  the  effect  of  sediment  rigidity  on  bottom  reflection  loss 
(RL)  from  typical  deep  sea  sediment  types  shows  that  sediment  shear  (S) 
wave  excitation  is  important  for  thin  sediment  layers  but  is  negligible 
for  thick  layers.  The  major  mechanism  for  S  wave  excitation  is  compress ional 
(P)  wave  conversion  at  the  sediment-substrate  interface:  Little  energy 
is  coupled  into  S  waves  at  the  water-sediment  interface. 

The  mechanism  for  S  wave  excitation  provides  a  means  for  quantifying 
the  categories  of  thick  and  thin  layers.  It  is  the  amplitude  of  the  P 
wave  at  the  substrate  interface  that  matters.  If  the  P  wave  amplitude 
at  the  substrate  interface  is  significant,  the  layer  is  thick;  if  it  is 

ill,  the  layer  is  thin.  In  general  thr  re  is  a  grazing  angle  separating 
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a  low  angle  region  in  which  the  P  wave  turning  point  is  well  above  the 
substrate  from  a  high  angle  region  in  which  the  P  wave  strikes  the  substrate 
interface  or  has  a  turning  point  near  the  interface.  For  angles  less  than 
0q  the  layer  is  thick  and  S  wave  effects  are  negligible.  For  angles  greater 
than  0^  the  layer  is  thin  and  S  wave  effects  produce  significant  increases 
in  RL.  For  physically  thick  layers  0q  is  large  and  S  wave  effects  are 
negligible  for  all  grazing  angles,  except  for  a  shift  in  the  peak  structure 
between  the  critical  angles  for  substrate  S  and  P  wave  propagation.  For 
physically  thin  layers  0^  can  be  0°  and  sediment  rigidity  will  be  particular¬ 
ly  important  at  low  grazing  angles. 

For  thin  layers  the  dependence  of  RL  on  frequency  shows  the  existence 
of  high  and  low  frequency  regimes.  The  total  attenuation  of  the  S  wave 
traveling  one  way  through  the  sediment  provides  a  means  of  separating  these 
frequency  regimes.  For  total  S  wave  attenuation  greater  than  20  dB  an 
S  wave  generated  at  the  substrate  interface  is  essentially  absorbed  within 
the  sediment  without  striking  the  other  interface.  Since  we  assume  attenuation 
is  proportional  to  frequency,  this  defines  a  high  frequency  regime  in  which 
the  propagation  of  sediment  S  waves  can  be  neglected;  however,  their  excitation 
is  still  an  important  loss  mechanism.  At  low  frequencies  S  wave  attenuation 
is  less  than  20  dB  across  the  layer  and  the  additional  interference  effects 
due  to  the  propagation  of  S  waves  through  the  sediment  are  important. 

Very  large  peaks  in  RL  ( >20  dB)  are  possible. 
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TABLE  I.  Acoustic  parameters  of  the  hypothetical  turbidite  layer: 

c  is  the  compression.^  wave  speed,  v  is  the  shear  wave  speed, 

P  is  the  density,  k^  is  the  compress ional  wave  attenuation, 

k  is  the  shear  wave  attenuation, 
s 
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FIGURE  3 

REFLECTION  LOSS  versus  FREQUENCY  FOR  A  252  m 
THICK  HYPOTHETICAL  TURBIDITE  LAYER  AT  20° 


ARL:UT 
AS-79-1815-P 
PJV-  GA 
8- 17-  79 


SOLID  SEDIMENT 
FLUID  SEDIMENT 


10 


APPENDIX  C 


A  RAY  PATH  ANALYSIS  OF  SEDIMENT  SHEAR  WAVE  EFFECTS 
ON  BOTTOM  REFLECTION  LOSS 


by 


Paul  J.  Vidmar 


Paper  submitted  to  The  Journal  of  The  Acoustical  Society  of  America 
for  publication 


105 


■  ***** 


A  ray  path  analysis  of  sediment  shear  wave  effects  on  bottom  reflection 
loss 
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The  effect  of  sediment  shear  wave  propagation  on  bottom  reflection 

loss  (RL)  from  a  typical  deep  sea  sediment  layer  is  analyzed  by 

means  of  a  ray  path  decomposition  of  the  acoustic  field.  The 

amplitude  of  eacii  shear  (S)  and  compressional  (P)  ray  is  obtained 

from  the  interface  reflection  and  transmission  coefficients. 

These  coefficients  are  analyzed  by  means  of  an  expansion  in 

the  parameter  e  which  is  proportional  to  the  relatively  small 

S  wave  velocity  of  marine  sediments.  Ray  paths  with  amplitudes 
2 

up  to  order  e  are  identified  and  used  to  develop  a  qualitative 

ray  path  model  of  RL.  P  wave  conversion  at  the  sediment- 

substrate  interface  is  the  dominant  mechanism  for  sediment  S  wave 

2 

excitation.  To  order  e  ,  RL  is  the  result  of  the  interference 
of  three  waves  in  the  water  whose  amplitudes  are  proportional 
to  e°,  e\  and  z" .  Comparison  with  RL  obtained  from  a  com¬ 
putational  model  shows  Chat  the  ray  path  model  correctly  predicts 
the  dependence  of  RL  on  frequency,  S  wave  speed,  S  wave  attenua¬ 
tion,  and  P  wave  at  enuation. 

PACS  number:  43.20.Fn  43.20.Dk  43. 30. Dr  43. 30. Bp 
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INTRODUCTION 


I. 

Recent  work  on  the  effect  of  sediment  rigidity  on  bottom  reflection 

loss,  RL  ■  -20  l°g^g  (|R|)»  where  R  is  the  reflection  coefficient  of  the 

ocean  bottom,  shows  that  sediment  shear  (S)  wave  excitation  can  be  an 

important  loss  process  influencing  low  frequency  sound  propagation  in 

the  oceans.  For  a  typical  deep  sea  sediment  structure,  a  single  turbi- 

dlte  layer  overlying  a  basalt  substrate,  the  sediment  S  wave  excitation 

1  2 

is  negligible  at  low  grazing  angles  for  thick  sediment  layers.  ’ 

However,  for  thin  layers  sediment  S  wave  excitation  is  an  important,  and 
in  some  cases  dominant,  energy  loss  process.^  The  dependence  of  RL  on 
frequency  for  thin  layers  has  low  and  high  frequency  regimes  with 
strikingly  different  properties.^  At  low  frequencies  RL  increases  a  few 
decibels  relative  to  the  RL  for  a  fluid  sediment  and  also  has  very  large 
peaks  up  to  25  dB.  At  high  frequencies  the  peaks  disappear,  and  rather 
broad  oscillations  in  RL  occur  in  both  fluid  and  solid  sediment  cases 
with  the  solid  sediment  RL  having  an  almost  constant  increase  above  the 
fluid  sediment  RL.  Compressional  (P)  wave  conversion  at  the  sediment- 
substrate  interface  was  suggested  as  the  major  mechanism  for  S  wave 
excitation  in  the  sediment.^ 

In  this  paper  a  ray  path  analysis  of  P  and  S  wave  propagation  in 
a  solid  sediment  layer  is  developed  and  used  to  elucidate  the  qualitative 
properties  of  RL  described  above.  Basic  to  this  approach  is  an  under¬ 
standing  of  the  effect  of  sediment  rigidity  on  the  reflection  and  trans¬ 
mission  coefficients  at  the  water-sediment  and  sediment-substrate 
interfaces.  An  expansion  of  these  coefficients,  using  a  small  parameter 

e  proportional  to  the  small  S  wave  speed  in  marine  sediments,  is  used 
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to  identify  P  and  S  wave  ray  paths  with  amplitudes  up  to  second  order  in 

the  expansion  parameter.  The  ray  path  model,  based  on  the  P  and  S  wave 

amplitudes,  reveals  the  mechanism  for  sediment  S  wave  excitation  to  be 

primarily  P  wave  conversion  at  the  sediment-substrate  interface  and  shows 

that  RL  can  be  described  as  the  result  of  the  interference  of  three 

waves  representing  energy  returned  to  the  water  from  the  sediment.  The 

three  waves  have  magnitudes  proportional  to  the  expansion  parameter  to 

the  zero,  first,  and  second  powers.  The  utility  of  the  ray  model  is 

established  by  comparing  predictions  of  the  dependence  of  RL  on  sediment 

3 

parameters  with  results  obtained  from  a  numerical  model  of  RL  from 
solid  sediments. 

The  remainder  of  this  paper  is  organized  as  follows.  In  section  II 
the  expansion  theory  of  the  reflection  and  transmission  coefficients 
is  developed.  In  section  III  these  coefficients  are  used  to  develop 
the  qualitative  ray  path  model.  In  section  IV  the  bottom  reflection 
loss  generated  by  a  numerical  model  is  used  to  check  the  predictions 
and  applicability  of  the  ray  model.  Section  V  discusses  the  mechanism 
for  shear  wave  excitation  in  sediments.  The  final  section  summarizes 
our  results. 

II.  EXPANSION  THEORY  OF  REFLECTION  AND  TRANSMISSION  COEFFICIENTS 

The  following  conventions  and  definitions  will  be  used  in  the 

study  of  the  reflection  and  transmission  coefficients.  Subscripts  1, 

2,  and  3  identify  quantities  in  the  water,  sediment,  and  substrate, 

respectively.  For  simplicity  we  assume  a  horizontally  stratified  ocean 

bottom  with  homogeneous  media,  and  a  harmonic  time  dependence  with 
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*  n 


angular  frequency  w. 

Compress ional 

wave  speeds 

are  c^ ,  S  wave  speeds 

are  ,  and  densities 

are  p ^ .  Wave 

numbers  are 

defined  as 

r  2  2l1/2 

Y“'Y  V"/vr 

k*KjCOse, 

Frkl  • 

[V  Brk]  • 

where  6  la  the  grazing  angle  of  the  Incident  wave  in  the  water. 

Subscripts  p  and  s  will  identify  reflection  and  transmission  coeffi¬ 
cients  for  P  and  S  waves  respectively.  The  water-sediment  interface 
will  be  denoted  by  W-S  and  the  sediment-substrate  interface  by  S-S. 

The  basis  of  our  study  of  the  reflection  and  transmission 
coefficients  is  an  expansion  in  the  small  parameter  This  is 

an  excellent  expansion  parameter  at  the  W-S  interface  where  the  surficial 
shear  wave  velocities^  are  between  50  m/sec  and  200  m/sec  giving 
0.03<e<0.15.  The  expansion  parameter  is  not  as  good  at  the  S-S  interface 

because  of  the  larger  S  wave  velocity  due  to  positive  shear  speed 

4  4 

gradients  in  realistic  sediments.  For  a  typical  turbidite  layer  e 

varies  from  0.076  at  the  W-S  interface  to  0.185  at  36  m,  0.256  at  120  m, 

and  0.406  at  518  m.  Even  at  518  m  the  second  order  terms  contribute 

only  15%  corrections  compared  to  the  40%  corrections  due  to  first  order 

terms.  Only  the  lowest  order  corrections  in  e  will  be  retained  in 

our  expansion. 

The  z  axis  Increases  downward  from  the  water  into  the  sediment. 

The  displacement  u^  is  defined  in  terms  of  a  scalar  potential, 

ei^X  an<*  a  vector  potential,  Xj“yXj(z)  ut\  by 

J7®j+.£x  .2j »  where  y  is  a  unit  vector  along  the  y  axis.  The  depth 
dependencies  are  given  by 
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1 


ix.z  -lie  z 

■  H"  e  J  +  W*  e  J  ,  and 


j  j 


i 


16.  z  -16. z 

Y  »  v”  e  ^  +  V*’  e  J 

J  j  j  * 


where  U  and  V  are  the  amplitudes  o£  the  downward  propagating  waves 

•J  J 

and  and  V.  are  the  amplitudes  of  the  upgoing  plane  waves. 

J  J 

Beginning  at  the  W-S  interface  with  a  unity  amplitude  sound  wave 

Ik.  z  —  lie.  z  .  ..  . 

incident  from  the  water,  $^**e  +^pi2e  U  '  and,  following 

the  procedure  given  in  the  appendix,  we  obtain 

Rpl2  R12  2e  T12  T21  COS  9  * 


(1) 


T  *  j.  2e2T^°^  cos^O 

pi 2  i2  +  T12  R21  COS  9 


(2) 


ts12  -  2e2t[°2)  (x2/k)  cos20 


(3) 


where  R^-C^Pj-K^p^/U^+K^)  and  T^°^-2Kipi/(Kipj+KjPi)  are  the 

zero  order  reflection  and  transmission  coefficients  obtained  for  a  fluid 

sediment  with  the  F  wave  incident  from  medium  i.^ 

Equations  (1)  through  (3)  show  that,  as  far  as  the  incident  sound 

wave  is  concerned,  the  sediment  can  be  accurately  treated  as  a  fluid. 

Equations  (1)  and  (2)  show  that  sediment  rigidity  produces  only  second 
2 

order  e  corrections  to  the  fluid-fluid  interface  reflection  and  trans¬ 
mission  coefficients.  Equation  (3)  shows  that  the  excited  sediment 

2 

shear  wave  amplitude  is  of  order  t  and  thus  contains  little  of  the 
incident  energy. 

The  downward  propagating  P  and  S  waves  next  strike  the  S-S  interface. 
For  the  incident  P  wave  we  find 
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<4> 


T 


Rp23  "  Rp23  “  erpn  0086 


Tp23  “  Tp23  "  erpTp  cose 


R 


S23  "  _2erp  C080  * 


Ts23  “  Tb23  ~  erpT8  C°80 


For  an  Incident  S  wave  we  have 


(5) 

(6) 

(7) 


rp23  "  n(l-eraco80) 


(8) 


Cp23  =  Tp(l-er8cos0) 


1  -  2eT  cos0  , 

a  9 


(9) 

(10) 


and 


ts23  "  Ts(1_er8cO80) 


(ID 


(16) 


i 


T 

8 


T(o)k 

i+W  + 


r<°> 

P23 


and  the  zero  order  reflection  and  transmlsalon  coefficients  for  a 
P  wave,  with  horizontal  wave  number  k,  incident  from  a  fluid  sediment:**'** 


Rp23  “  {p3K2[(1~2a2)2  +  "  P2K3}/D  ’  (17) 

Tp23  ”  2p2<2(l“2a  )/d  ,  (18) 

lg23  m  ^p2K2a  k3  I  ^  >  (19) 

where  a-^/c^cosG  and 

D  *  p3<2  [(1_2a2)2  +  4a4<2e3 /k2]  +  P2IC3  *  (20) 


Equations  (4)  through  (11)  show  that  sediment  rigidity  can  be 

important  at  the  S-S  interface.  According  to  Eq.  (6)  the  P  wave  incident 

upon  the  S-S  interface  excites  a  first  order  (e^*)  S  wave  in  the  sediment. 

Equations  (4),  (5),  and  (7)  show  that  the  reflected  P  wave  and  P  and 

S  waves  transmitted  into  the  substrate  also  have  first  order  corrections 

to  their  amplitudes .  Recall  that  the  S  wave  excited  at  the  W-S  inter- 
2 

face  is  of  order  e  [Eq.  (3)].  Equations  (8)  through  (10)  show  that 
the  reflected  and  transmitted  waves  due  to  it  are  at  most  of  order  e2 
and  are  small  compared  to  the  larger,  order  e^-,  corrections  due  to  P 
wave  conversion  at  the  S-S  interface.  The  presence  of  gradients  will 
Increase  the  relative  effectiveness  of  the  S-S  interface  by  increasing 
the  value  of  e. 

Equations  (8)  through  (11)  show  an  apparent  inconsistency.  The 

incident  S  wave  is  perfectly  reflected  to  order  e°  and  yet  sets  up  other 
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reflected  and  transmitted  waves  whose  amplitudes  are  also  of  order  e°. 

A  consideration  of  the  total  energy  flux  resolves  this  inconsistency. 
Recalling  that  the  energy  flux  is  proportional  to  the  inverse  of  the 

wave  speed, the  incident  and  reflected  S  waves  in  the  sediment  have 

1 

an  energy  flux  proportional  to  c  while  the  energy  flux  of  the  other 
waves  is  to  lowest  order  proportional  to  e°.  The  S  waves,  then,  carry 
one  order  higher  energy  flux  than  the  other  waves  even  though  all  the 
amplitudes  are  of  the  same  order.  To  order  e  \  the  incident  and 
reflected  S  waves  conserve  energy  flux.  To  the  next  order,  c°,  the  part 
of  the  S  wave  energy  that  is  not  perfectly  reflected  (e*  corrections  to 
the  amplitude)  provides  some  order  e°  energy  which  appears  in  the  lowest 
order  energy  flux  of  the  other  waves. 

We  next  consider  incident  S  and  P  waves  in  the  sediment  striking 
the  W-S  Interface.  The  computed  reflection  and  transmission  coeffi¬ 
cients  are 


p21 

-  +  2e2Tj°)T2°)cos20 

9 

(21) 

p21 

-  T*J*  +  Re2!^*^00*2® 

9 

(22) 

ls21 

-  2e2T^)(K1/k)cos20 

(23) 

for  an  incident  P  wave,  and 

rp21  '  -2c(1-R21))CO‘9  •  <24) 


',21  ■  2cTnc°*6 

rs21  “  -1 


(25) 

(26) 
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1 


for  an  incident  S  wave.  The  lowest  order  corrections  to  Eq.  (26)  are 
3 

of  order  c  . 

Equations  (21)  through  (23)  give  the  same  qualitative  results  as 

Eqs.  (1)  through  (3).  Both  sets  of  equations  show  that  the  corrections 

2 

due  to  sediment  rigidity  are  of  order  e  at  the  W-S  interface  when  a 

P  wave  is  incident,  l.e.,  the  sediment  can  be  accurately  treated  as  a 

fluid  at  the  W-S  interface  whenever  a  P  wave  is  incident.  Recalling 

that  the  incident  S  wave  is  at  most  order  e[Eqs.  (3),  (6),  (10)],  the 

corrections  to  P  wave  amplitudes  due  to  S  wave  conversion  are  at  most 
2 

of  order  e  [Eqs.  (24),  (23)].  This  leads  to  the  conclusion  that,  as 
far  as  P  waves  are  concerned,  the  sediment  at  the  W-S  interface  can  be 
accurately  treated  as  a  fluid. 

III.  RAY  PATH  MODEL  AND  PREDICTIONS 

The  transmission  and  reflection  coefficients  derived  above  will 

now  be  used  to  develop  a  ray  path  model  of  the  qualitative  effects  of 

sediment  rigidity  in  a  single  homogeneous  sediment  layer.  The  reflection 

and  transmission  coefficients  of  the  previous  section  will  be  used  to 

determine  the  wave  amplitudes  associated  with  the  ray  paths  at  each 

Interface.  Subscripts  U  or  L  will  identify  the  values  of  the  expansion 

parameters  at  the  W-S  or  S-S  interfaces,  respectively.  The  corrections 

due  to  sediment  rigidity  at  each  interface  can  then  be  identified  and 

the  qualitative  effect  of  the  shear  speed  gradient  treated  by  setting 

2 

eT >e  .  To  allow  for  e  >e  ,  corrections  of  order  eT  will  be  retained 
L  u  L  u  L 

2 

while  those  of  the  order  eu  will  be  neglected  compared  to  one. 

The  ray  path  model  is  shown  schematically  in  Figs.  1,  2,  and  3. 

The  solid  lines  are  ray  paths  for  P  waves  and  the  dashed  lines  are  those 
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for  S  waves.  Each  ray  path  is  labeled  with  the  amplitude  of  the  wave 
associated  with  that  path.  The  positions  at  which  the  rays  strike  the 
interfaces  are  labeled  a,  b,  c,  and  d  for  Identification  purposes. 

Figure  1  shows  the  ray  paths  excited  when  the  S  wave  attenuation  is  high 
enough  for  the  S  wave  to  be  totally  absorbed  before  crossing  the  layer. 
Figure  2  shows  the  additional  ray  paths  resulting  from  the  first  reflec¬ 
tion  of  the  sediment  S  waves,  and  Fig.  3  shows  those  due  to  the  second 
reflection.  Since  the  S  wave  speed  in  marine  sediments  is  so  small 
compared  to  the  sound  speed  in  water,  the  S  wave  ray  paths  are  almost 
normal  to  the  plane  of  stratification. 

In  Fig.  1  the  S  wave  attenuation  is  high  enough  for  an  S  wave  to 
be  completely  absorbed  in  traveling  one  way  through  the  sediment.  At 
position  a,  the  incoming  wave  in  the  water  reflects  from  the  W-S  inter¬ 
face  with  negligible  corrections  [Eq.  (1)]  due  to  sediment  rigidity. 

The  corrections  to  the  transmitted  P  wave  are  also  negligible  [Eq.  (2)]. 

2 

A  small  amplitude  S  wave  of  order  eu  is  also  excited  [Eq.  (3)].  At 

position  b  the  transmitted  P  wave  strikes  the  S-S  interface  where  an 
.  2 

order  eL+eL  ®  wave  is  excited  [Eq.  (6)].  The  reflected  P  wave  and  the 

2 

substrate  P  and  S  waves  suffer  order  e^+e^  corrections  [Eqs.  (4),  (5), 

and  (7)].  The  reflected  P  wave  carries  its  corrected  amplitude  back 

through  the  sediment  to  the  W-S  interface  at  position  c.  Here  a  small 
2 

order  S  wave  is  excited  [Eq.  (23)].  The  transmitted  wave  in  the 
2 

water  has  ej+e^  corrections  [Eq.  (22)].  The  reflected  P  wave  carries 

2 

corrections  of  the  order  eL+eL  (Eq.  (21)]  back  down  through  the  sediment 
where  it  again  interacts  with  the  S-S  interface  at  position  d,  and  the 
process  repeats. 
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Figure  2  shows  Che  additional  rays  set  up  when  the  S  wave  attenuation 


is  reduced  to  allow  the  sediment  S  waves  to  strike  one  interface  before 

2 

being  absorbed.  The  order  S  wave  ray  originating  at  the  W-S  interface 

(position  a  in  Fig.  1)  is  totally  reflected  from  the  substrate  at 

2 

position  a  [Eq.  (10)]  and  retains  its  eu  magnitude  (corrections  of 
2  2 

order  e  e.  are  neglected).  Order  e  P  and  S  waves  are  also  excited  in 
u  L  u 

the  substrate  at  position  a  [Eqs.  (9)  and  (11)]  and  a  reflected  P  wave 
2 

of  order  e  in  the  sediment  [Eq.  (8)].  At  position  b  the  P  wave  from 

2 

position  a  transmits  its  magnitude  into  the  water  and  also  into  the 
reflected  P  wave  at  the  W-S  interface  [Eqs.  (21)  and  (22)].  The  reflected 
S  wave  is  of  order  and  is  neglected  [Eq.  (23)].  At  position  c  the 

downward  propagating  P  wave  strikes  the  S-S  interface  where  it  sets  up 

2  2 
order  eu  waves  in  the  substrate  and  an  order  eu  reflected  P  wave 

[Eqs.  (4),  (5),  and  (7)].  No  significant  S  wave  is  excited  in  the 

sediment  [Eq.  (6)].  The  upgoing  P  wave  then  continues  to  position  d 

where  the  process  begun  at  position  b  repeats. 

Figure  2  also  contains  ray  paths  coming  from  the  S  wave  originating 

at  the  substrate  interface  (position  b  in  Fig.  1).  At  position  b  this 

S  wave  is  totally  reflected  from  the  W-S  interface  [Eq.  (26)]  and  retains 

its  order  eT+e^  magnitude.  A  reflected  P  wave  of  order  eueL  excite<^ 

[Eq.  (24)]  along  with  a  wave  of  amplitude  water  I^q.  (25)]. 

The  reflected  P  wave  carries  its  magnitude  through  a  reflection  at  the 

S-S  interface  at  position  c  [Eq.  (4)]  where  it  also  excites  P  and  S 

waves  of  order  e  et  in  the  substrate  [Eqs.  (5)  and  (7)].  No  significant 

sediment  S  wave  is  set  up  [Eq.  (6)].  The  upward  traveling  P  wave  is 

again  reflected  from  the  W-S  interface  at  position  d.  The  reflected 
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sediment  P  wave  and  the  transmitted  wave  in  the  water  carry  an  order 
eueL  magnitude  [Eqs.  (21)  and  (22)].  No  significant  S  wave  is  excited 
(Eq.  (23)].  The  order  e  eT  P  wave  moving  downward  from  the  position  d 
then  continues  the  process. 


Figure  3  shows  the  situation  for  an  attenuation  small  enough  for  the 

sediment  S  waves  to  reflect  a  second  time  before  being  absorbed.  The 

S  waves  originating  at  the  W-S  interface  (positions  a  and  c  in  Fig.  1) 

produce  negligible  effects  when  they  strike  the  W-S  interface  at  positions 

2 

a  and  c.  The  reflected  S  waves  retain  their  order  ey  magnitudes  and  no 

significant  P  waves  are  excited  [Eqs.  (24),  (25),  and  (26)].  The  S 

wave  originating  at  the  S-S  interface  (position  b  in  Fig.  1),  however, 

does  produce  additional  corrections  to  the  energy  returned  to  the  water. 

On  striking  the  S-S  interface  at  position  b  it  reflects  with  its  order 
2 

£j+£l  magnitude  [Eq.  (10)].  It  also  excites  P  and  S  waves  in  the  sub- 

2 

strate  and  a  reflected  P  wave  in  the  sediment,  all  of  order  e  +e  [Eqs.  (8), 

L  L 

(9) ,  and  (11) ] .  The  P  wave  travels  through  the  sediment  to  the  water 

2 

interface  where  it  reflects  at  position  c,  retaining  its  order  eT+eT 

Lt  L 
2 

magnitude,  and  is  transmitted  into  the  water  with  magnitude  eT+£. 

L  L 

[Eqs.  (21)  and  (22)].  It  does  not  excite  a  significant  S  wave 

[Eq.  (23)].  At  position  d,  the  downward  propagating  P  wave  reflects 

2 

from  the  substrate  interface  retaining  its  eT+eT  magnitude,  and  sets  up 

it  L 

2 

order  waves  in  the  substrate  [Eqs.  (4),  (5),  and  (7)].  The 

reflected  P  wave  then  continues  the  process. 

2 

Figure  3  also  contains  an  S  wave  of  order  excited  in  the  sediment 
at  position  d  [Eq.  (6)].  This  S  wave  produces  no  further  effects  since 
it  is  absorbed  before  reaching  the  W-S  interface,  as  can  be  seen  by 


recalling  that  the  magnitude  of  this  S  wave  is  proportional  to  the 
amplitude  of  the  P  wave  which  struck  the  S-S  interface.  The  magnitude  of 
the  P  wave  is  in  turn  proportional  to  the  magnitude  of  the  S  wave  which 
originated  at  the  S-S  interface  and  has  traveled  up  and  back  through 
the  sediment.  The  maximum  amplitude  of  the  newly  excited  S  wave  then 
corresponds  to  the  amplitude  of  the  twice  reflected  shear  wave.  By 


hypothesis  (in  this  example)  an  S  wave  is  absorbed  before  its  third 
reflection.  Hence,  the  new  S  wave  will  be  absorbed  before  it  reaches 
the  W-S  interface. 

For  still  lower  S  wave  attenuation  additional  reflections  from  the 
interfaces  are  important.  New  ray  paths  similar  to  those  of  Figs.  2  and 


3  are  generated.  The  lower  S  wave  attenuation  means  that  the  ray 

2 


paths  of  a  particular  type  (for  example  those  of  order  simi^ar 

to  those  of  Fig.  3)  have  comparable  amplitude  but  differ  in  phase 
because  of  the  additional  path  lengths  traveled  by  the  S  waves. 

For  a  homogeneous  layer  (e^e^e)  Figs.  1  through  3  show  that  there 
are  three  types  of  waves  emerging  from  the  W-S  interface.  Figure  1 


shows  the  largest  of  these.  The  directly  reflected  ray  and  the  series 


of  rays  carrying  a  magnitude  of  order  1+e+e  combine  to  form  a  single 


wave,  the  A.  wave.  For  a  fluid  sediment  (t=0)  the  A  wave  still  exists 
o  o 


and  in  fact  produces  the  usual  fluid  sediment  RL.  For  e^O  the  Aq  wave 


contains  a  correction  to  its  magnitude  of  order  e+e  due  to  the 
S  wave  excited  at  the  substrate  interface  by  the  P  wave.  The  phase  of 
the  A  wave  does  not,  however,  depend  on  S  wave  propagation.  The  second 


type  of  wave,  the  A^  wave,  is  formed  by  the  series  of  rays  of  order 


e+e  shown  in  Fig.  3.  It  is  the  most  important  correction  to  the  A 
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wave,  ics  magnitude  is  proportional  to  e  ,  to  lowest  order,  while  its 

phase  contains  the  two-way  S  wave  phase  ^  The  third  type  of 

2 

wave,  the  wave,  is  formed  by  the  series  of  rays  of  order  e  in 

2 

Fig.  2.  Since  its  magnitude  is  of  order  e  it  produces  only  small 

corrections  compared  to  those  due  to  the  wave.  The  phase  of  the  A2 

wave  contains  the  one-way  S  wave  phase, 

The  ray  path  model  predicts  the  dependence  of  the  magnitudes  of 

the  A^  and  A^  wave  on  shear  wave  attenuation.  Figures  2  and  3  show 

that  the  magnitude  of  these  waves  depends  on  the  attenuation  suffered 

by  the  S  wave  in  traveling  through  the  sediment.  The  A^  wave  (Fig.  3) 

depends  on  the  two-way  attenuation  while  the  wave  (Fig.  2)  depends 

only  on  the  one-way  attenuation.  As  S  wave  attenuation  increases,  the 

A^  wave  will  decrease  in  magnitude  more  rapidly  than  the  Az  wave. 

For  low  S  wave  attenuation,  powerful  interference  effects  are 

possible  between  the  many  ray  paths  of  a  particular  type  (for  example, 

2 

additional  e  paths)  returning  energy  to  the  water.  These  interference 
effects  will  modify  the  dependence  of  the  A^  and  A^  waves  on  S  wave 
attenuation.  The  magnitude  of  Aq  is  also  affected  since  essentially  all 
of  the  energy  converted  from  P  to  S  must  return  to  the  water. 

Finally,  the  ray  path  model  shows  the  importance  of  S  wave 
excitation  by  P  wave  conversion  at  the  substrate  interface.  All  correc¬ 
tions  of  order  e  occur  because  of  the  initial  interaction  of  the  P  wave 
with  the  substrate  shown  in  Fig.  1.  The  order  e  corrections  in 
the  Aq  wave  are  clearly  due  to  this  interaction  (Fig.  1).  The  A^  wave 
is  excited  by  S  waves  set  up  by  the  initial  P  wave  interaction  (Fig.  3). 
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Part  of  the  second  order  corrections,  the  e  e  contribution  in  Fig.  2, 

U  Li 

is  also  due  to  this  initial  P  wave.  The  remaining  part  of  A2»  that 
2 

labeled  cu  in  Fig.  2,  while  basically  set  up  by  P  wave  conversion  at  the 
water  interface,  also  depends  on  the  P  wave  traveling  through  the  sediment 
layer.  The  same  P  wave  conditions  which  would  make  the  wave 

unimportant  would  also  affect  this  contribution  to  A2  in  the  same  manner. 
Thus,  it  is  crucial  for  the  P  wave  amplitude  at  the  S-S  interface  to  be 
significant  if  sediment  rigidity  is  going  to  have  an  important  effect 
on  RL. 

IV.  COMPARISON  OF  RAY  PATH  MODEL  PREDICTIONS  WITH  COMPUTED  BOTTOM 
REFLECTION  LOSS 

In  this  section  a  computational  model,  described  in  detail 
3 

elsewhere,  will  be  used  to  test  some  of  the  predictions  of  the  ray  path 
model.  The  dependence  of  R  on  frequency  that  is  produced  by  the  com¬ 
putational  model  will  be  used  to  extract  the  magnitudes  of  the  Aq,  A^, 
and  A2  waves.  Homogeneous  layers  will  be  treated  at  frequencies  near 

14  Hz  and  a  grazing  angle  of  20°.  The  acoustic  parameters  used  in  the 

3 

numerical  calculations  are:  p^  ®  1.053  g/cm  ,  c^  =  1540  ra/sec, 

3  3 

P2  =  1.27  g/cm  ,  c2  =  1525  ra/sec,  v2  =  70  m/sec,  =  2.8  g/cm  , 

c^  =  5700  m/sec,  and  v^  =  2900  ra/sec.  The  P  wave  attenuations  are: 

k  »  =  0.115  dB/ra/kHz,  k  _  =  0.152  dB/m/kHz.  The  S  wave  attenuations 
p2  p3 

are:  kg2  =  10.0  dB/m/kHz,  kg^  =  0.90  dB/m/kHz.  The  layer  thickness 
is  H  =  40  m.  Individual  parameters  will  be  varied  to  illustrate  the 
features  of  the  ray  path  model.  The  range  of  the  parameters  does  not 
necessarily  represent  values  attained  in  marine  sediments. 
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The  different  phases  of  the  A^,  A^,  and  A ^  waves  sake  it  possible 
to  extract  their  amplitudes  from  the  interference  structure  seen  in  the 
dependence  of  RL  on  frequency.  The  ray  path  model  shows  that  the 
magnitude  of  the  reflection  coefficient  can  be  written  approximately  as 

|R|  *  Aq  +  A1  cos y £  +  A2  cosCt^+u)  ,  (27) 

where  Aq,  A^,  and  A2  are  the  amplitudes  of  the  Aq  (order  1+e) , 

1  2 

A^  (order  e  )  and  A2  (order  c  )  waves,  and  v  is  a  frequency  independent 

phase  difference  between  the  and  A2  waves.  This  form  for  ) R ]  assumes 

that  0^b,.jH/v2  and  $2*u2H /v2  contain  the  major  frequency  dependence  of 

the  phases.  This  is  true  in  our  case  since  the  frequency  range  will  be 

limited  so  that  the  phase  shifts  due  to  P  wave  propagation  and  the 

interface  continuity  conditions  are  small  compared  to  those  due  to  S 

wave  propagation.  The  major  interference  pattern  in  j R j .  will  be  due  to  a 

superposition  of  the  Aq  and  A1  waves.  The  smaller  magnitude  A2  wave  will 

add  a  modulation  to  this  structure. 

A  typical  reflection  coefficient  is  shown  in  Fig.  4.  Four  of  the 

maxima  and  minima  are  labeled  by  letters.  The  basic  oscillatory  structure 

is  due  to  the  interference  of  the  A  and  A,  waves.  The  maxima  (or  minima) 

o  1 

recur  when  ^  changes  by  2t.  The  difference  in  height  of  peaks  b  and  d 
clearly  shows  the  effect  of  the  modulation  produced  by  A2>  The  unknown 
amplitudes  and  the  relative  phase  in  Eq.  (27)  are  obtained  by  setting: 
?1*;2=0  at  peait  a*  c2a"  at  valley  b;  at  peak  c; 

;^*=3_/2,  ?2=^"r  at  valley  d.  This  generates  a  set  of  four  equations: 


R  =  A  +  A.  +  An  cos; 
a  o  1  2 

» 

(28) 

R^  *  A  -  A^  -  A2  sin; 

• 

(29) 
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R  »  A  +  A-  -  A-  cos$ 
c  o  1  2 


(30) 


Rd  “  Ao  ~  A1  +  A2  sin*  *  (31) 

where  R^,  R^,  R^,  and  R^  are  Che  measured  values  of  |r|  at  a,  b,  c, 
and  d.  These  four  equations  can  be  easily  solved  for  the  four  unknowns. 

The  magnitudes  of  Aq,  ,  and  A£  obtained  from  Eqs.  (28)  through 
(31)  are  not  sufficiently  accurate  for  the  purposes  of  this  section. 

The  reason  for  this  inaccuracy  is  an  inconsistency  in  Eqs.  (28)— (31) ; 
the  frequency  dependence  of  Aq,  A^,  and  A2  has  not  been  included,  while 
R^,  R^,  Rc,  and  Rd  were  measured  at  different  frequencies.  In  practice  it 
is  easier  to  correct  the  procedure  for  measuring  [ R |  than  to  modify 
Eqs.  (28)-(31).  A  more  accurate  procedure  for  measuring  |R|  is  also 
illustrated  in  Fig.  4.  Peaks  corresponding  to  similar  interference 
maxima  (or  minima)  are  joined  by  straight  lines:  a  a',  b  b',  and  d  d'. 

The  values  of  R^,  R^,  and  are  then  estimated  at  the  frequency  of 
maxima  c  by  obtaining  their  magnitudes  from  the  intersection  of  a 
vertical  line  through  maxima  c  with  lines  aa' ,  bb',  and  dd'.  The 
resulting  values  for  Ra>  R^,  Rc>  and  now  have  the  major  frequency 
dependence  removed  and  can  be  used  in  Eqs.  (28)-(31)  to  obtain  more 
accurate  values  of  the  magnitudes  of  Aq,  A^,  and  A2.  This  more  accurate 
procedure  was  used  to  obtain  the  results  of  this  section. 

The  first  parameter  we  will  investigate  is  the  sediment  shear  speed 
v 2 •  Figure  5  shows  the  computed  RL  for  3  different  shear  speeds: 

150  m/sec,  70  m/sec,  and  a  fluid  case,  V2=0  m/sec.  The  results  are  in 
qualitative  agreement  with  the  ray  model  predictions.  The  RL  generally 
increases  with  sediment  shear  speed.  The  rapid  oscillatory  structure, 
due  to  the  interference  of  the  Aq  and  A^  waves,  increases  in  magnitude 
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and  spreads  out  in  frequency  as  V2  increases.  The  nodulation  due  to  the 
small  A£  wave  is  clearly  evident  in  the  V2  •  150  m/sec  curve.  A 
similar  oscillatory  structure  was  observed  in  the  low  frequency  regime 
of  dependence  of  RL  on  frequency  computed  for  a  hypothetical  turbidite 
layer.1 

Figure  6  shows  the  dependence  of  the  frequency  interval  between 
adjacent  maxima  of  RL,  Af,  on  sediment  shear  speed.  The  line  is  the 
ray  path  model  prediction,  Af*V2/2H.  The  agreement  is  excellent. 

Figure  7  shows  the  dependence  of  the  amplitudes  of  the  A^  and  A2 

waves  on  sediment  shear  speed.  The  third  set  of  points  shows  the  change 

of  A  due  to  sediment  rigidity,  i.e.,  AA  (e)**A  (0)-A  (c).  For  reference, 
o  000 

v2  ”  150  m/sec  corresponds  to  £“0.10.  In  agreement  with  the  ray  path 
model  predictions,  A^  and  are  proportional  to  V2  at  small  shear 
speeds.  The  Increase  in  AAq  is  consistent  with  an  increase  in  RL, 
which  the  ray  path  model  predicts  will  be  proportional  to  e  to  lowest 
order.  All  three  sets  of  points  show  the  effect  of  higher  order 
corrections  at  larger  e.  The  magnitude  of  A2  is  much  smaller  than  that 
of  A^,  in  agreement  with  the  ray  path  model  predictions.  The  lines  in 
Fig.  7  are  drawn  to  illustrate  the  proportionality  to  e;  the  slopes  do 
not  reflect  theoretical  predictions. 

Figures  2  and  3,  and  the  previous  discussion  of  the  ray  path  model, 
show  that  the  magnitude  of  the  and  A2  waves  have  different  depend¬ 
encies  on  the  S  wave  attenuation,  i.e.,  A^«exp(-2k^gH)  and  kf&W (-ik^H) , 
where  kig  is  the  imaginary  part  of  the  S  wave  wave  number.  Figure  1  shows 
that  the  magnitude  of  the  Aq  wave  is  independent  of  S  wave  attenuation. 

The  dependence  of  Aq,  A^,  and  A2  on  the  sediment  shear  wave  attenuation 
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Is  shown  In  Fig.  8.  The  abscissa  in  Fig.  8  is  k^,  the  S  wave  attenuation 
in  dB/m/kHz.  Since  we  assume  a  linear  dependence  of  k^g  on  frequency, 
kg  is  related  to  k^g  by  k^g*kgf /8686.  The  values  of  kg  greater  than  10 
attenuate  the  S  wave  by  about  6  dB  in  a  distance  H.  This  effectively 
eliminates  any  waves  other  than  those  shown  in  Figs.  1  through  3  and 
permits  the  ray  model  to  predict  a  clean  exponential  dependence  on  kg 
for  the  magnitudes  of  the  A^  and  waves.  The  lines  in  Fig.  8  are 
"eyeball"  fits  to  the  data  obtained  from  the  computed  reflection  coeffi¬ 
cients.  The  slopes  are  0.0285  for  the  A^  wave  and  0.0569  for  the  A^ 
wave,  and  they  differ  by  a  factor  of  2  as  predicted  by  the  ray  theory  and 
agree  very  well  with  the  predicted  slopes  of  0.028  and  0.056.  The 
change  in  Aq  due  to  S  waves  AAq  is  seen  to  be  independent  of  kfi  as 
expected. 

The  ray  path  model  also  predicts  that  AAq,  A^,  and  depend  on 

the  attenuation  of  the  P  wave  in  the  sediment.  Figures  1  through  3  can 

be  used  to  deduce  this  dependence  for  large  P  wave  attenuations. 

Figure  1  shows  that  AAq  is  carried  by  a  P  wave  traveling  down  and  back 

through  the  sediment,  AAQ«exp(-2kipH) .  Figure  2  shows  that  A^exp^k^H) 
2 

also.  The  ey  part  of  clearly  has  this  dependence.  The  eu£L  part  has 
the  same  dependence  since  it  comes  from  the  initial  P  wave  striking 
the  substrate  Interface  at  position  b  in  Fig.  1.  Figure  3  combined  with 
Fig.  1  shows  that  A^*exp(-2k^pH) .  These  relations  apply  at  higher  P  wave 
attenuations  where  only  a  single  ray  path  returning  to  the  water  has  a 
significant  amplitude. 

Figure  9  shows  the  dependence  of  AA  ,  A. ,  and  A.  on  P  wave 

o  1  L 

attenuation.  The  abscissa  in  Fig.  9  is  k  ,  the  P  wave  attenuation,  in 

P 
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dB/m/kHz.  The  velues  of  were  conputed  from  the  actuel  Imaginary  part 

of  the  wave  number.  Including  corrections  due  to  grazing  angle.  For 

this  study  the  sediment  thickness  is  H  ■  80  m.  The  values  of  k  above 

P 

5  attenuate  the  ?  wave  by  more  than  6  dB  in  a  distance  H.  This  effectively 
eliminates  all  but  the  first  ray  returning  from  the  sediment  to  the  water 
in  Figs.  1  through  3.  The  ray  path  model  then  predicts  the  exponential 
dependence  given  above.  The  lines  are  again  "eyeball"  fits  to  the  cal¬ 
culated  magnitudes.  The  slopes  of  the  lines  are  0.059  for  A^,  0.123  for 
A^,  and  0.133  for  AAq.  The  A^  and  Aq  slopes  are  approximately  equal  and 
are  nearly  double  the  slope  of  the  A^  wave  in  agreement  with  the  ray 
model.  The  values  of  the  slopes  also  agree  well  with  the  predicted 
values  of  0.056  and  0.112. 

The  ray  path  model  strictly  predicts  no  important  S  wave  effects 

if  the  sediment  P  wave  speed  is  larger  than  the  sound  speed  in  water. 

The  incoming  wave  in  the  water  is  specularly  reflected  with  only  order 
2 

t  corrections.  The  computational  model,  based  on  wave  theory,  allows 
the  penetration  of  the  P  wave  into  the  sediment  where  it  becomes 
evanescent  and  decays  exponentially.  Generalizing  the  ray  path  model 
prediction  for  the  case  of  P  wave  attenuation  to  the  case  of  the 
evanescent  P  wave,  one  obtains  the  same  predictions  with  k^,  now  given  by 

kip  "  2irf  [<c2cos0/cx)2  -  l]  /c2 

Since  this  exponential  decay  constant  is  proportional  to  frequency,  it 

can  be  written  in  terms  of  a  frequency  Independent  decay  constant, 

k',  in  units  of  dB/m/kHz,  i.e.,  k.  -k'f/8686. 

P  IP  P 
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Figure  10  shows  the  dependence  of  AAQ,  Aj ,  and  A 2  on  k\  For  this 

study  H  “  100  m.  For  the  P  wave  amplitude  decays  by  about  11  dB  in 

a  distance  H.  The  parameter  actually  varied  in  this  case  was  the  P  wave 

velocity  in  the  sediment.  The  lines  are  "eyeball"  fits  to  the  calculated 

magnitudes.  The  measured  slopes  are:  0.162  for  AA  ,  0.162  for  A,,  and 

o  1 

0.072  for  The  slope  of  the  A^  line  is  approximately  half  that  of 

the  slopes  of  the  AAq  and  A^  lines  in  agreement  with  the  ray  path  model 
predictions.  The  magnitude  of  the  slopes  also  agrees  well  with  the 
predicted  values  of  0.070  and  0.140. 

The  qualitative  effects  of  a  gradient  in  shear  speed  can  be  deduced 
from  Figs.  1  through  3.  The  gradient  causes  t.  to  be  larger  than  c  ,  and 
Figs.  1  and  3  show  that  the  larger  increases  the  effect  of  sediment 
rigidity  on  the  Aq  and  A^  waves.  Since  these  two  waves  dominate  the 
structure  of  the  reflection  loss,  the  increase  in  RL  due  to  sediment 
rigidity  will  be  larger  in  layers  with  a  shear  speed  gradient  than  in  a 
homogeneous  layer.  Figure  2  shows  that  a  larger  e  also  increases  the 
cueL  part  of  the  A^  wave,  making  it  possible  for  the  shear  speed  at  the 
W-S  Interface  to  slightly  influence  RL;  this  would  occur  at  low  fre¬ 
quencies  where  the  interference  of  the  A-  wave  with  the  A  and  A,  waves 

i  ox 

is  important. 

The  use  of  the  actual  reflection  and  transmission  coefficients  given 
above  could  be  used  to  develop  a  quantitative  ray  path  model.  The  pro¬ 
cedure  would  be  analogous  to  the  local  plane  wave  approach  used  recently^ 
to  develop  a  quantitative  ray  path  model  of  fluid  sediments. 
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VI.  MECHANISM  FOR  SHEAR  WAVE  EXCITATION  IN  DEEP  SEA  SEDIMENTS 


Coupreaslonal  wave  conversion  at  the  sediment-substrate  Interface 

has  been  proposed  as  the  major  mechanism  for  shear  wave  excitation  in  a 

typical  deep  sea  sediment.^  The  expansion  theory  of  the  reflection  and 

transmission  coefficients  and  the  ray  path  model  developed  here  support 

this  conclusion.  The  generality  of  the  assumptions  of  the  expansion 

theory  and  ray  path  analysis  shows  that  this  mechanism  in  fact  applies 

to  a  large  class  of  marine  sediments.  The  only  restrictions  are  that  the 

sediment  be  a  single,  inhomogeneous  layer  and  that  gradient-driven  P-S 

coupling  be  negligible.  For  layered  rather  than  continuous  sediment 

structures,  additional  conversion  at  Interfaces  within  the  sediment  may 

also  be  important.  For  sand  type  sediments,  important  in  shallow  water 

continental  shelf  environments,  the  large  near-surface  gradients^’  may 

make  the  gradient-driven  P-S  coupling  near  the  water-sediment  interface 

9 

an  important  additional  mechanism  at  low  frequencies.  For  typical  deep 

sea  sediments,  the  gradient-driven  P-S  conversion  is  probably  negligible 
3 

above  about  3  Hz. 

For  typical  deep  sea  sediment  types,  the  mechanism  for  sediment  S 

wave  excitation  to  lowest  order  in  c  can  be  deduced  from  the  combined 

action  of  the  reflection  and  transmission  coefficients.  As  far  as  P 

waves  are  concerned  the  sediment  is  a  fluid  at  the  water-sediment  interface, 

and  negligible  energy  is  transferred  into  S  waves.  An  order  t  S  wave  is 

excited  in  the  sediment  by  the  P  wave  striking  the  sediment-substrate 

interface.  This  sediment  S  wave  has  no  further  effect;  it  is  essentially 

perfectly  reflected  from  the  water-sediment  and  sediment-substrate  inter- 

2 

faces  and  produces,  at  most,  order  e  corrections  to  the  P  wave  returned 
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to  the  water  column.  The  P  wave  reflected  from  the  sediment-substrate 
Interface  does,  however,  have  an  Important  role  to  play;  Its  amplitude 
picks  up  an  order  e  correction  at  the  sediment-substrate  interface.  This 
correction  travels  through  the  water-sediment  interface  and  affects  the 
energy  returned  to  the  water. 

The  following  picture  emerges.  Sediment  rigidity  affects  bottom 
reflection  loss  by  means  of  energy  transferred  to  sediment  shear  waves 
at  the  sediment-substrate  interface.  Sediment  compressional  waves  play 
a  major  role;  they  excite  the  shear  wave  at  the  sediment-substrate  inter¬ 
face  and  are  responsible  for  carrying  the  influence  of  the  shear  wave 
back  through  the  sediment  and  into  the  water.  The  major  effect  of  sediment 
rigidity  is  then  to  provide  an  energy  sink  at  the  sediment-substrate 
interface. 

The  empirical  classification  of  sediments  into  thick  and  thin 
categories  based  on  the  effect  of  sediment  shear  waves  can  now  be  under¬ 
stood  and  quantified.  For  a  thick  sediment  layer,  the  P  wave  has  its 
turning  point  well  above  the  sediment-substrate  interface  and  it  does 

not  significantly  interact  with  the  sediment-substrate  interface.  Since 
2 

only  order  e  corrections  are  possible  due  to  the  refracted  P  waves  at 
the  water-sediment  interface,  sediment  rigidity  has  a  negligible  effect 
on  bottom  reflection  loss.  In  this  case  the  sediment  can  be  accurately 
treated  as  a  fluid.  For  a  thin  sediment  layer,  the  P  wave  interacts 
significantly  with  the  sediment-substrate  interface  where  sediment  S  waves 
are  excited,  and  the  sediment  must  be  treated  as  a  solid.  Whether  a 
sediment  layer  is  thick  or  thin  is  then  translated  into  the  more  easily 
quantified  question  of  whether  the  P  wave  interacts  with  the  sediment- 


substrate  interface  (thick)  or  not  (thin).  The  answer  to  this  question 
depends  on  frequency  (through  P  wave  attenuation)  and  grazing  angle 
(because  of  P  wave  evanescence) . 

An  additional  important  conclusion  of  this  paper  is  that  sediment 

rigidity  can  be  neglected  at  the  water-sediment  interface.  The  ray  model 

predicts  that  the  lowest  order  effects  depending  on  the  S  wave  velocity 

2 

at  the  water-sediment  interface  are  contained  in  the  c  and  e  e.  waves 

u  u  L 

seen  in  Fig.  2.  Since  e'vO.l  these  are  small  corrections  compared  to  the 
corrections  seen  in  Figs.  1  and  3  due  to  P  wave  conversion  at  the 
sediment-substrate  interface. 

The  mechanism  for  sediment  S  wave  excitation  also  leads  to  a 
generalization  of  the  "hidden  depths"  concept^  of  fluid  sediment  theory. 
According  to  this  concept  there  is  a  depth  below  which  the  subbottom 
structure  has  a  negligible  effect  on  the  acoustic  field  in  the  water. 

This  depth  is  essentially  a  few  wavelengths  below  the  turning  point  of 
the  P  wave  in  the  subbottom.  Below  this  turning  point  the  P  wave  decays 
exponentially.  The  acoustic  energy  interacting  with  the  structure  below 
the  hidden  depth  is  a  very  small  fraction  of  the  total  energy.  Hence, 
the  structure  below  the  hidden  depth  is  acoustically  negligible.  Since 
sediment  S  waves  originate  principally  by  means  of  the  sediment  P  wave 
interacting  with  the  sediment-substrate  interface,  the  hidden  depth 
criterion  also  applies  to  sediment  S  wave  effects.  If  the  sediment- 
substrate  interface  is  in  the  hidden  depths,  the  P  wave  amplitude  is 
negligible  at  the  sediment-substrate  interface  and  sediment  S  wave  effects 
are  negligible.  For  solid  sediments,  the  structure  in  the  hidden  depths 
is  still  negligible  and,  in  addition,  if  the  substrate  interface  is  in 
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the  hidden  depths,  the  effect  of  sediment  S  waves  is  also  negligible  and 
the  sediment  can  be  accurately  treated  as  a  fluid. 

Figure  11  illustrates  the  relative  insensitivity  of  R  to  the 
sediment  S  wave  velocity  at  the  water-sediment  interface  (V2U)  compared 
to  that  at  the  sediment-substrate  interface  (v2L).  Figures  11(a)  and 
11(b)  use  the  computational  model  to  follow  the  development  of  |r| 
starting  with  a  homogeneous  layer  having  v^  *  100  m/sec  and  ending  with 
V2  *  300  m/sec.  To  eliminate  interference  effects  a  high  frequency 
(S  wave  heavily  attenuated)  was  chosen.  In  Fig.  11(a)  V2u  is  constant 
while  v^L  increases.  A  constant  gradient  is  used  which  has  the  value  of 
5  sec  ^  when  v^L  ■  300  m/sec.  Three  grazing  angles  are  shown.  The 
| R |  decreases  almost  linearly  with  V2L  at  each  grazing  angle.  The 
decrease  is  significant,  30%  at  8*30°  and  50%  at  8=10°.  Figure  11(b)  shows 
the  additional  change  in  | R  j  which  occurs  when  V2L  is  held  constant  at 
300  m/sec,  V2^  is  increased  to  fora  a  homogeneous  layer,  and  with 
V2  *  300  m/sec.  As  expected,  the  factor  of  3  increase  in  V2L  now  pro¬ 
duces  only  small  changes  in  |r|  compared  to  those  produced  by  changing 
V2L*  This  result  reinforces  the  conclusion  that  RL  is  relatively 
insensitive  to  the  surficial  S  wave  velocity.  In  fact,  fairly  accurate 
values  of  |r|  can  be  obtained  at  high  frequencies  by  ignoring  the 
S  wave  velocity  gradient  and  assuming  a  homogeneous  layer  with  V2 
given  by  the  S  wave  velocity  at  the  sediment-substrate  interface.  For 
the  same  gradients  used  in  Fig.  11(a),  the  error  in  |r|  is  less  than 
12%  over  the  entire  range  of  gradients  and  for  grazing  angles  less  than 
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50°.  This  leads  Co  the  conclusion  that  v.  is  the  important  S  wave 

ZL 

parameter  at  high  frequencies. 


VI .  SUMMARY 

A  theoretical  treatment  of  the  reflection  and  transmission 
coefficients  at  the  water  and  substrate  interfaces,  using  an  expansion 
in  the  small  sediment  S  wave  velocity,  provides  a  basis  for  developing 
a  ray  model  of  the  effect  of  sediment  rigidity  on  bottom  reflection 
loss.  The  basic  mechanism  for  S  wave  excitation  emerges  as  P  wave 
conversion  at  the  sediment-substrate  interface.  The  dependence  of  RL 
on  S  wave  velocity  and  attenuation,  and  P  wave  attenuation,  are  pre¬ 
dicted  by  the  ray  model  and  verified  in  detail  by  results  obtained 
from  a  computational  model.  The  bottom  reflection  loss  is  found  to  be 

due  to  the  interference  of  three  waves:  the  A  wave,  similar  to  that 

o 

generated  by  a  fluid  sediment,  the  A^  wave,  whose  amplitude  depends  on 
the  S  wave  velocity  to  the  first  power,  and  the  A^  wave,  whose  amplitude 
depends  on  the  S  wave  velocity  to  the  second  power. 
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APPENDIX 


The  procedure  for  calculating  the  reflection  and  transmission 
coefficients  to  any  order  in  e  is  illustrated  in  this  appendix.  The 
results  obtained  for  a  P  wave  incident  upon  the  S-S  interface  will  be 
used  as  an  example.  (The  notation  of  section  II  is  used.) 

The  continuity  of  the  horizontal  and  normal  components  of 
displacement  and  the  normal  and  shear  stress  can  be  written  in  matrix 
form  as 


Mx  =  y 


(A-l) 


(+  -  +  -  \T 

W2*  W3*  ^2*  ^3/  contains  the  reflected  and  transmitted 
potential  amplitude;  Xc,[“P2(l-^a2j»-K2,^P2a2IC2  / lc»— k/B2J  »  the  stress" 
displacement  vector  due  to  the  incident  unity  amplitude  P  wave;  and 


-P3(1-2^) 

2p2a2S2/k 

2p3a3B: 

1  -‘2 

-<3 

k 

-k 

M  * 

1  v/fj* 

2p3a32<3/k 

°2(1-2*2) 

-P3(l- 

V/e2 

-k/62 

1 

e3/B 

(A-2) 


in  which  aj  =  (v_. /c^)cos0 .  The  superscript  T  denotes  the  transpose  of  the 
row  vector.  The  quantities  8^  anc*  a2  in  —  anci  X  contain  the  expansion 
parameter  e=V2/c^.  The  rows  of  M  correspond  to  the  normal  stress,  normal 
displacement,  shear  stress,  and  horizontal  displacement,  respectively. 

The  vectors  x  and  and  the  Matrix  M  are  now  expressed  as  a  series 
in  e  as 


x  =  x  +  ex  +  e  x„  +...+ 


(A- 3) 
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(A-4) 


+V 


1  m  io  +  ^  +  E  Z2  +•  *  •  + 

=  "  =o  +  %.  +  e2==2  +‘  •  *+  »  (A-5) 

Substituting  Eqs.  (A-3)-(A-5)  into  Eq.  (A-l)  and  equating 
quantities  of  the  same  order  in  e  yields  a  sequence  of  equations; 

^  ,  (A-6) 

M^i  -  v i  -  Xq  .  (A"7) 

=0^2  “  *2  ~  =2^o  “  =l-i  ‘  (A_8) 

These  equations  can  be  solved  in  sequence:  Eq.  (A-6)  yields  x^, 

Eq.  (A-7)  and  x  yield  x_  »  Eq.  (A-8)  along  with  x  and  jc  yields  x«* 

U  1  O  i.  A 

The  quantities  in  M  and  that  depend  on  e  are  a2**ecos0,  and 
3  3 

k/82=ecos0-e  cos  0/2+...+  .  The  quantities  in  Eqs.  (A-6)-(A-8)  are  then: 


-P3(l+a*) 


3  3 
0 


V" 


0  -2p, 


M, 


COS0 


°\ 


0  0 

0  0 

1  -1 


0 

0 


0 

0 


(A-9) 


(A- 10) 


(A-ll) 


o  e3/k 
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yx  -  -cose<0,  0,  0,  I)1 


y2  "  2P2COS^e^1’  °»  K2^k’ 


The  solution  of  Eqs.  (A-6)-(A-8)  using  Eqs.  (A-9)-(A-14)  is 
straightforward  and  yields  Eqs.  <4)-(7)  in  section  II.  The  process 
can  be  continued  to  any  order  in  e. 
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COMPR ESSIOH AL  WAVE 


- SHEAR  WAVE 


RAY  PATHS  EXCITED  WHEN  THE  SEDIMENT  SHEAR  WAVE 
IS  ABSORBED  BEFORE  REFLECTION 
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COMPRESSION*!.  WAVE 


ADDITIONAL  RAY  PATHS  DUE  TO  THE  FIRST  REFLECTION 
OF  THE  SEDIMENT  SHEAR  WAVES 
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REFLECTION  LOSS  -  dB 


FIGURE  5 

REFLECTION  LOSS  rtrsu*  FREQUENCY  FOR 
A  40  m  THICK  HOMOGENEOUS  LAYER  WITH 
SEDIMENT  SHEAR  SPEED  AS  A  PARAMETER 
FOR  CURVE  (a)  v2  =  ®  (b)  v2  =  30  m/»«e; 

(c)  v2  =»  ISO  m/*»c 
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FREQUENCY  SEPARATION  -  Hi 


FIGURE  6 

THE  DEPENDENCE  OF  THE  FREQUENCY  INTERVAL 
BETWEEN  THE  ADJACENT  MAXIMA  OF  RL 
ON  SEDIMENT  SHEAR  SPEED 
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50  100 

SHEAR  SPEED  -  m/*«c 


FIGURE  7 

THE  DEPENDENCE  OF  THE  MAGNITUDES  OF  aA 
A,,  AND  A2  ON  SEDIMENT  SHEAR  SPEED 


r 


5 

COMPRESSIONAL  WAVE  ATTENUATION  -  dB/fa/kH* 


FIGURE  9 

THE  DEPENDENCE  OF  THE  MAGNITUDES  OF  A  An, 
A,,  AND  A2  ON  THE  ATTENUATION  COEFFICIENT 
OF  THE  SEDIMENT  COMPRESSIONAL  WAVE,  kp 
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DECAY  CONSTANT  -  dB/m/kHz 


FIGURE  10 

THE  DEPENDENCE  OF  THE  MAGNITUDES  OF  aA0, 
A,,  AND  A2  ON  THE  DECAY  CONSTANT  OF  THE 
EVANESCENT  SEDIMENT  COMPRESSIONAL  WAVE,  kj, 
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FIGURE  11 

THE  DEPENDENCE  OF  THE  MAGNITUDE  OF  THE  REFLECTION 
COEFFICIENT  ON  THE  SEDIMENT  SHEAR  WAVE  SPEED  AT 
INTERFACES  FOR  GRAZING  ANGLES  OF  10°  (CIRCLES), 

20°  (TRIANGLES),  AND  30°  (SQUARES) 

(o)  v 2y  CONSTANT,  (b)  v2L  CONSTANT 
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